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ABSTRACT 


Langmuir  circulations  in  a  coastal  environment  at  an  inner  shelf  site  with  strong 
tidal  foreing  were  examined  using  a  4-month  observation  of  high-resolution  veloeity 
profiles  during  the  ONR-sponsored  CBLAST-Low  air-sea  interaetion  experiment. 
Because  of  their  potential  eontributions  to  mixing  proeesses  in  the  surface  boundary 
layer,  Langmuir  eireulations  are  important  for  naval  operations  such  as  ASW,  MCM  and 
Amphibious  operations  in  littoral  waters.  A  deteetion  algorithm  for  upwelling  / 
downwelling  veloeities  due  to  Langmuir  eireulations  was  developed.  Analyses  focused 
on  long-feteh,  unstratified  water  eolumn  eonditions  with  loeally  developed  seas.  The 
observations  showed  that  strong  Langmuir  eireulations  greatly  redueed  water  eolumn 
shear  arising  from  the  tidally  foreed  bottom  boundary  layer,  while  there  was  no 
signifieant  effeet  from  the  strong  tidal  boundary  layer  on  the  formation  of  Langmuir  eells. 
Long  feteh  conditions  had  the  strongest  vertieal  veloeities  with  Langmuir  number  0. 1-0.4, 
while  short  fetch  conditions  have  weaker  vertical  velocities  with  Langmuir  number 
greater  than  0.5.  Cell  vertieal  velocities  were  correlated  with  the  water  frietion  veloeity, 
but  penetration  depths  did  not  depend  on  foreing  mechanisms.  Spacing  and  penetration 
depths  were  eorrelated  for  inereasing  wind  speed  and  wave  heights  whereas  this 
eorrelation  was  not  seen  for  eonstant  high  wind  speed  and  wave  heights. 
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I.  INTRODUCTION 


The  physical  processes  controlling  atmosphere  ocean  momentum  and  scalar 
fluxes  in  the  presence  of  waves  remain  an  important  research  topic  in  Physical 
Oceanography  and  Meteorology.  Meteorologists  are  trying  to  understand  the  effects  of 
ocean  waves  in  high  and  low  wind  conditions  by  examining  the  temperature,  turbulence 
and  humidity  profiles  in  the  Marine  Boundary  Layer  (MBL)  while  the  oceanographers 
are  interested  in  measuring  surface  fluxes  in  Ocean  Boundary  Layers  (OBL). 

Mixing  processes  in  the  wind  forced  ocean  surface  boundary  layer  (OSBL)  are 
very  different  from  the  wall  boundary  layers  due  to  the  action  of  surface  gravity  waves 
that  modify  the  momentum,  heat  and  gas  exchange  between  the  MBL  and  OSBL.  The 
main  physical  processes  that  govern  the  mixing  in  OSBL  are  convective  forcing  due  to 
night  time  cooling,  shear  stress  instability  generated  by  wind  stress  acting  through  waves 
and  wave  breaking,  and  Langmuir  circulations  resulting  from  the  interactions  between 
surface  gravity  waves  and  the  wind  driven  surface  currents. 


A.  MOTIVATION  OF  THIS  STUDY 

Langmuir  circulations  in  a  coastal  environment  are  investigated  in  this  study  to 
understand  the  role  of  coherent  structures  seen  below  wind-forced  waves  in  mixing 
across  the  ocean  mixed  layer.  In  inner  shelf  coastal  environments,  the  surface  layer 
frequently  interacts  with  a  tidally  forced  bottom  boundary  layer  in  unstratified  conditions. 
It  is  important  for  modern,  high-resolution  regional  numerical  models  to  correctly 
parameterize  these  multiple  turbulent  processes  in  littoral  regions. 

Langmuir  circulations  (Figure  1.1)  are  counter  rotating  vortices  with  convergence 
zones  above  the  downwelling  regions  and  divergence  zones  above  the  upwelling  regions 
at  the  sea  surface.  They  exist  in  the  surface  boundary  layer  due  to  the  interactions 
between  strong  wind  driven  surface  currents  and  the  Stokes  drifts  of  surface  gravity 
waves.  Under  unstratified  conditions,  continental  shelf  circulations  can  penetrate  down  to 
the  ocean  bottom  and  have  dominant  effects  on  the  turbulent  kinetic  energy  budget 
(TKE)  across  the  water  column. 
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Figure  1.1 


Wind 


The  mixing  processes  in  a  coastal  region.  The  eddies  due  to  surface 
mixing  processes  in  surface  boundary  layer  and  the  eddies  due  to  strong 
tidal  flow  in  bottom  boundary  layer  are  separated  by  the  seasonal 
stratification  (summer  time)  and  merge  with  unstratified  ocean  (winter 
time).  MLD  is  the  mixing  layer  depth. 


Langmuir  circulations  in  a  coastal  environment  with  strong  tidal  forcing  were 
examined  using  a  4-month  observation  of  high-resolution  velocity  during  CBLAST 
(Coupled  Boundary  Layers  Air  Sea  Transfer)  at  a  coastal  observatory  south  of  Martha’s 
Vineyard  during  the  fall  of  2003.  The  study  focuses  on  conditions  with  long  fetch  waves, 
an  unstratified  ocean  water  column  and  locally  well-developed  waves.  An  algorithm  to 
detect  downwelling  /  upwelling  regions  of  Langmuir  circulations  was  created  for  the 
vertically  oriented  beam  of  five  beam  BADCP  (Broadband  Acoustic  Doppler  Current 
Profiler)  high-resolution  velocity  profiles.  A  review  of  ocean  TKE  budget  and  Langmuir 
circulations  overview  is  given  in  Chapter  11.  Method  of  data  collection  with  five  beam 
BADCP  and  description  of  CBLAST  experiment  site  is  introduced  in  Chapter  Ill. 
Methods  used  for  data  processing  to  detect  Langmuir  circulations  and  their  features  are 
explained  in  Chapter  IV.  Analysis,  including  definition  of  the  supporting  conditions  of 
Langmuir  circulations,  Langmuir  number  calculation  and  relationship  with  the  detected 
downwelling  /  upwelling  velocities,  the  correlations  between  the  downwelling  / 
upwelling  and  the  water  friction  velocity,  horizontal  and  vertical  length  scales  is 
discussed  in  Chapter  V.  Conclusions  are  presented  in  Chapter  VI. 
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II.  TURBULENCE  FLOW  AND  LANGMUIR  CIRCULATIONS 


A.  TURBULENCE  FLOW 

Turbulent  flow  exists  in  the  upper  oeean  due  to  the  mixing  proeesses  foreed  by 
eonveetion  and  free  eonveetion.  In  this  study,  foreed  eonveetion  due  to  surfaee  wind 
stress  in  the  presenee  of  surfaee  gravity  waves  will  primarily  be  eonsidered.  The 
important  parameters  to  define  the  turbulent  flow  are  the  Riehardson  number  and  the 
Reynolds  numbers.  In  this  ehapter,  the  turbulent  kinetie  energy  equations  for  stratified 
flow  are  summarized  to  provide  an  overview  of  turbulent  proeesses  found  in  the  eoastal 
oeean. 


The  gradient  shear  stability  Riehardson  number  is  a  measure  of  the  stability  of  the 
oeean  flow,  and  it  is  the  ratio  of  statie  stability  and  dynamie  stability. 


Ri  =  - 


shear  I 


(2.1) 


where  N  is  Brunt-Vaisala  or  buoyaney  frequeney. 


N=-^x 


S  SPe 


p0  Sz 


(2.2) 


dp 

The  g  is  gravitational  aeeeleration, - -  is  the  stratifieation  measured  by  the  vertieal 

dz 

density  gradient  and  z  is  the  depth,  whieh  is  positive  upward.  Negative  vertieal  density 
gradients  represent  a  hydrostatieally  unstable  oeean,  whieh  lead  to  eonveetive  instability, 
while  a  positive  vertieal  density  gradient  represent  a  hydrostatieally  stable  oeean. 
Sustained  hydrostatieally  unstable  eonditions  arise  from  free  eonveetion  forees  that  are 
typieally  due  to  thermal  eooling  of  the  oeean  surfaee. 


The  shear  term  in  the  denominator  of  Equation  2.1  is  another  faetor  responsible 
for  turbulent  flow  and  represents  the  foreed  eonveetion.  Foreed  eonveetion  in  the  eoastal 
oeean  is  due  to  wind  foreing,  internal  wave  interaetions,  frontal  and  meso-seale 
horizontal  shear,  and  the  eurrent  interaetions  with  the  seabed.  Shear  in  the  water  eolumn 
ereates  horizontal  and  vertieal  vortiees.  If  the  Riehardson  number  is  less  than  %, 
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turbulent  flow  is  more  likely  to  exist.  The  higher  the  shear  term,  the  smaller  the 
Riehardson  number  is,  and  the  stronger  the  turbulent  flow. 


z 


V 


Figure  2.1  x,  y  and  z  coordinate  system  with  vorticity  circles. 


For  the  coordinate  system  defined  in  Figure  2.1,  u,  v  and  w  are  the  x,  y,  z  velocity 
components  of  the  velocity  field  respectively.  The  vorticitiy  equations  for  x,  y  and  z-axis 
are  given  in  the  Equations  2.3  through  2.5  respectively. 


dv  dw 
CO  = - 1 - 

dz  dy 


(2.3) 


du  dw 
dz  dx 


(2.4) 


du  dv 
dy  dx 


(2.5) 


All  the  terms  on  the  right  hand  side  of  the  Equations  2.3  through  2.5  represent  fluid  shear 
components.  The  z-vorticity  Equation  2.5  doesn’t  have  any  vertical  component  of  the 
flow,  which  means  there  is  no  vertical  mixing  in  this  equation.  The  vertical  mixing 
process  exists  only  in  the  x-vorticity  and  y-vorticity  equations. 

Reynolds  number  is  another  estimation  of  strength  of  turbulent  flow,  defined  as 
the  ratio  of  inertial  forces  (advection  terms)  and  viscous  forces  in  the  equations  of  motion 
(Equation  2.6) 

Re  =  -^  (2.6) 

V 
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Where  U  is  the  veloeity  seale,  L  is  the  horizontal  length  seale  and  o  is  the 
viseosity.  If  Reynolds  numbers  are  mueh  larger  than  one  or  greater  than  2000,  turbulent 
flow  is  likely  to  exist  due  to  negligible  moleeular  frietion  effeets.  This  eondition  is  true 
not  only  for  turbulent  flow,  but  also  for  every  kind  of  flow  in  the  oeean.  The  frietion  term 
in  the  equations  of  motion  (momentum  equations)  depends  upon  the  fluid  eharaeteristies 
and  for  most  seales  of  motion  in  the  oeean  ean  be  negleeted,  although  all  motions  are 
dissipated  due  to  the  fluid  viseosity  at  small  (0.1-10  mm)  seales. 


Turbulenee  needs  a  eontinuing  souree  of  energy  to  be  sustained.  The  basie 
balanee  of  turbulent  kinetie  energy  is  given  by 


d  (  TKE  ^  _  WORK  RA TE  VICOUS  DISSIPA TION 
dt  t  VOLUME )  ^  VOLUME  VOLUME 


(2.7) 


where  TKE  is  turbulent  kinetie  energy,  work  rate  is  eombination  of  the  work  done  by  free 
eonveetion  (thermally  driven)  and  foreed  eonveetion  (wind  foreed),  and  the  viscous  term 
is  the  dissipation  rate  of  the  turbulence,  which  is  commonly  known  as  s  . 

The  driving  properties  of  motion  in  the  ocean  can  be  generalized  into  the  three 
processes  of  momentum,  heat  content  and  salt  concentration.  The  motion  driven  by  these 
processes  have  mean  and  fluctuating  parts.  If  the  driving  property  is  called  F,  the 
conservation  laws  lead  to  Equation  2.8  and  2.9.  Equation  2.9  is  known  as  the  continuity 
equation. 


dE  BE  BE  BE  BE 

- = - VU - hV - \-W - =  /Cj. 

dt  Bt  Bx  By  Bz 


( B^E  B^E  B^e] 
^  Bx^  By^  Bz^  j 


+  Ep 


(2.8) 


Bu  Bv  Bw  _ 
—  +  —  + — =0 
Bx  By  Bz 


(2.9) 


Rp  represents  the  nonconservative  processes  in  (2.8).  The  mean  part  of  Equation  2.8  can 
be  combined  with  Equation  2.9  to  give  Equation  2.10, 


BE  -BE  -BE  —BE  BuE  BvE  BwE 

- h  U - h  V - h  W - 1 - 1 - 1 - =  Krj 

Bt  Bx  By  Bz  Bx  By  Bz 


B^E  B^E  B^E 
- - - 

Bx^  By^  Bz^ 


+  Rp  (2.10) 
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The  next  step  is  to  derive  the  variance  equation  for  F .  Variance  is  defined  as 


F'^  =  F' F' ,  this  variance  is  a  measure  of  the  turbulent  generation  of  fluctuations  inF  . 
By  subtracting  (2.10)  from  (2.8),  multiplying  by  2F' ,  and  taking  the  mean  of  each  term 
while  assuming  horizontal  homogeneity,  a  generalized  variance  budget  for  atmospheric 
and  oceanic  boundary  layers  results  as. 


d_ 

dt 


-zi^r^dF  d 

F  w - 

dz  dz 


F  ^ 


w 


dx 


+ 


dF 
V  ^ 


dz 


+  FR^ 


(2.11) 


The  first  term  is  the  unsteady  rate  of  increase/decrease  and  usually  called  “storage 
term”,  the  second  term  is  “gradient  production”  and  is  always  positive  because  of  the 
second  law  of  thermodynamics,  the  third  term  is  the  turbulent  diffusion  of  variance 
through  the  water  column,  the  fourth  term  is  the  rate  of  dissipation  variance  and  can  only 
be  negative,  the  last  term  is  collection  of  the  nonconservative  processes  if  they  exist. 


If  F  represents  a  horizontal  velocity  component  u  ,  the  turbulent  kinetic  energy 
equation  for  this  component  is. 


— 

d 

.2  \ 

u 

.  .  du 

, 

u 

=  -u  w  — 
dz 

dz 

w 

T , 

V  ; 

V 

/ 

u 


(2.12) 


f  a. .'A 


du 
ydxj 


f  a. . '  A 


+ 


du 

j 


+ 


^  du^ 

^dz  j 


(2.13) 


However,  the  processes  and  mathematical  terms  represented  by  R  have  not  been 
defined  yet.  The  u  -momentum  equation  will  be  used  to  make  this  definition. 


du  du  du  du 

- \-u - hv - \-w —  =  u 

dt  dx  dy  dz 


f  ^^2  ^2  ^2  ^ 

0  u  o  u  cu 


^'dz^ ) 


+  2Q^v-2Q  W-  — ^ 
p  dx 


(2.14) 


and  by  comparing  (2.14)  with  (2.8),  the  nonconservative  term  is. 
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1  dp 

R,  =2Qv-2Qw - ^ 

"  "  ^  p  dx 


(2.15) 


Since  all  the  terms  are  linear  in  Equation  2.15,  averaging  rules  can  be  used  to  get  the 
mean  version  of  R.  , 


Ru=2Qv-2Qw--^ 
^  p  dx 


(2.16) 


The  complete  mean  u  -momentum  equation  is, 


du  -du  -du  —du  duu  dvu  dwu 

- h  U - h  V - h  W - 1 - 1 - 1 - 

dt  dx  dy  dz  dx  dy  dz 


V 


^  d^u  d^u  d^u^ 


-  +  — — ^-  +  - 


+  2Q^v-2Q^,w-—^ 
p  dx 


(2.17) 


The  last  step  is  to  define  the  nonconservative  part,  F  R^  ,  of  the  variance  budget 
(2. 1 1)  by  using  (2.16), 


^ — r  .(  ,  .1  dp  ^ 

uR  u  =  u\  2Q  V  - 2Qw - — 


1  .  dp 


=  20.  uv  -20  uw - u  — 

p  dx 


V  P  SXy 

Substituting  (2.18)  in  (2.12)  the  final  version  of  the  u  -TKE  equation  is  obtained. 


(2.18) 


dt 


_ 

/  ,2 

u 

<  <du  d 

, 

u 

— 

=  -u  w - 

w 

— 

2  , 

dz  dz 

2  j 

V  7 

\  J 

p  du 


-s  2  - 20  uw  +  20  uv - 

“  p  dx 


(2.19) 


The  corresponding  Reynolds  and  TKE  equations  for  v  resemble  their  counterpart 
Equations  2.17  and  2.19,  except  for  the  Coriolis  terms. 


dv  -  dv  -dv  —dv  duv  dvv  dwv 

- h  U - h  V - h  W - 1 - 1 - 1 - 

dt  dx  dy  dz  dx  dy  dz 


u 


(  ^2  ^2  ^2 

c;  V  c  V  0  V 

-  + - ;r  +  - 


+  20 

p  dy 


(2.20) 
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(2.21) 


V 

dv  d 

f 

'  V 

— 

=  -vw - 

w 

1  - 

2  , 

dz  dz 

2  , 

V  ; 

V  / 

V 


-IQ^uv 


p  dv 
P  ^ 


The  vertical  balance  of  forces  in  flux  form  and  multiplied  through  by  density  is 
defined  in  Equation  2.22,  before  deriving  the  vertical  Reynolds  equation  and  vertical 
TKE  equation. 


P 


dw  duw  dvw  dww 

- 1 - 1 - 1 - 

dt  dx  dy  dz 


-  ^  dp 

=  2pau-- - pg  +  pv 

dz 


f  a2 


d  w  d  w  d  w 


-  +  - 


-  +  - 


(2.22) 


The  hydrostatic  approximation  “  0  = - pg  ”,  which  is  acceptable  for  large-scale 

dz 

motions,  neglects  the  vertical  accelerations  due  to  the  fluctuations  in  density  and  pressure 
and  will  exclude  the  turbulent  motion  and  small-scale  internal  waves.  It  is  necessary  to 
keep  the  vertical  forces  that  cause  the  vertical  accelerations  in  Equation  2.22  in  order  to 
describe  turbulent  motion.  The  “Boussinesq  Approximation”  can  be  used  to  accomplish 

this.  The  total  pressure,  p  ,  is  broken  into  two  parts,  hydrostatic  pressure  P  and  dynamic 
pressure  p  according  to  the  Boussinesq  approximation.  Hydrostatic  pressure  is  the 
pressure  that  would  be  caused  by  representative  or  average  density  p^  and  if  the  fluid 
were  motionless  with  no  horizontal  forces,  is  defined  as, 

dP 

^  =  -Pog  (2.23) 

dz 

Substituting  (2.23)  into  total  pressure  and  gravity  terms  in  (2.21), 


d  p 
dz 


Pg  = 


d{^P  +  p^ 


dz 


-Pg  =  g{Po-P)-^ 
dz 


(2.24) 


Buoyancy  force”  or  “reduced  gravity  force”  is  defined  as, 

g{Po-  Pit 


g  =b  = 

Substituting  (2.25)  in  (2.24), 


Po 


(2.25) 
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(2.26) 


^  r  ^  j 

- Pg  =  Pob 

oz 


dp 

dz 


Finally  substituting  (2.26)  in  (2.21),  getting  an  average  of  eaeh  term  and  dividing  eaeh 
term  by  in  (2.21),  the  vertieal  Reynolds  equation  is  obtained  as. 


dw 

dt 


duw  dvw  dww  duw  dvw  dww 


dx  dy 


2D.  u  - 


1  dp 
Po  dz 


+  b  +  o 


dz  dx 
d^w  d'^w 


dy 

d^w 


dz 


dx^  a/  dz^ 


(2.27) 


Buoyaney  foree  deseribed  in  (2.25)  has  mean  and  fluetuating  parts,  b  =  b  +  b  and 
b  =-p  g/  p^.  Fluetuations  in  buoyaney  foree  eause  positive  or  upward  and  negative  or 

downward  aeeelerations  attributable  to  the  negative  density  fluetuations  and  positive 
density  fluetuations  respeetively. 


w-TKE  equation  is  obtained  by  subtraeting  (2.27)  from  (2.22),  multiplying  by  w 
and  averaging  the  resultant  equation. 


d 

w 

5 

5 

w  p 

dt 

dz 

w 

I2J 

dz 

_  Po  _ 

P  dw 

+  bw  +2Duw  - s  2 

p^  dz  w 


(2.28) 


The  two  differenees  between  vertieal  TKE  Equation  2.28  and  the  horizontal  TKE 
equations  (2.21),  (2.19)  are  that  there  is  no  gradient  produetion  term  and  there  is  an  extra 

buoyaney  flux, bw  ,  'm  vertieal  TKE  equation.  There  is  also  no  gradient  produetion  term 
beeause  there  is  no  shear  in  w  if  mean  flow  is  horizontal. 


Summarizing  what  has  been  done  so  far,  three  eomponents  of  the  TKE  and  the 
Reynolds  equations  are  derived  assuming  horizontal  homogeneity.  The  u ,  v  and  w 
eomponents  of  TKE  are  rewritten  in  Equations  2.29,  2.30  and  2.31  respeetively  in  order 
to  explain  the  physieal  meanings  of  eaeh  term. 


— 

u 

<  •  du  d 

, 

u 

— 

=  -u  w - 

w 

I  - 

2  , 

dz  dz 

2  , 

V  7 

V  / 

- 2iluw  +2Q  wv 

p  dx  ^ 


(2.29) 
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_e 

dt 


— 

/  .2 

V 

dv  d 

, 

V 

— 

=  -v  w - 

w 

— 

2  , 

dz  dz 

2  j 

V  ; 

\  J 

P  “ 

-s  2 - 2Lluv 

”  p 


(2.30) 


_e 

dt 


( 

■  ^ 

w 

=  bw - 

, 

I  w 

p 

— 

w 

i - 1 - 

dz 

12 

Po) 

p  dw 

-s  2 - \-2ilu  w 

w  dz 


(2.31) 


The  physical  processes  identified  by  each  term  are  the  most  necessary  part  of 
understanding  the  turbulent  flow.  The  first  terms  are  the  time  rate  of  change  of  the  TKE 
or  “storage  terms”.  If  there  is  a  statistical  steady  state,  these  terms  are  equal  to  zero.  Due 
to  the  usual  rapid  time  adjustment  of  turbulence  to  changes  in  forcing,  these  terms  are 
often  assumed  to  be  negligible  compared  with  the  production  and  dissipation  terms 
described  below. 


The  second  terms  are  called  “Production  Terms”  which  represents  the  production 
processes  of  turbulence.  The  second  terms  in  (2.29)  and  (2.30)  are  the  shear  production 
terms  which  are  caused  by  forced  convection  and  the  Reynolds  stresses  working  against 
the  mean  shear.  The  second  term,  buoyancy  flux  term,  in  (2.31)  is  called  buoyant 

production  only  if  it  is  positive.  The  buoyancy  flux,  b'w  ,  will  be  positive  when  the  mean 
vertical  density  gradient  is  positive,  and  the  mean  vertical  buoyancy  gradient  is  therefore 
negative.  This  is  the  hydrostatically  unstable  case  and  free  convection  occurs  where 
consecutive  buoyant  parcels  of  fluid  rise  and  denser  particles  sink  due  to  the  gravity.  On 
the  other  hand,  whenever  the  mean  density  profile  is  stable  or  negative,  the  buoyancy  flux 
term  in  (2.31)  will  be  negative.  In  this  case,  the  vertical  mixing  due  to  mechanical  or 
forced  convection  will  be  damped  by  the  buoyancy  flux  term,  so  the  buoyancy  flux  term 
will  act  as  an  obstacle  to  the  vertical  TKE,  and  will  result  in  a  loss  of  the  vertical  TKE 
instead  of  a  gain  in  it.  This  is  frequently  the  case  for  the  oceanic  mixed  layer  in  the 
afternoon  due  to  solar  heating  of  the  ocean  surface  while,  the  marine  atmospheric 
boundary  layer  is  less  likely  to  become  strongly  stable  unless  warm  air  is  horizontally 
advected  over  cold  water  and  creates  the  inversion  layer.  Over  the  land,  the  atmospheric 
boundary  layer  is  likely  to  be  stable  due  to  the  radiation  inversion  at  night,  after  the 
ground  has  given  up  heat  by  long  wave  radiation  to  sky. 
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The  third  terms  represents  the  vertieal  transport  of  turbulenee  by  itself.  If  w 
earries  the  TKE,  adveetive  turbulent  transport  oeeurs.  The  vertieal  TKE  equation  has  an 
additional  kind  of  transport,  whieh  is  ealled  as  pressure  transport.  The  pressure  transport 
oeeurs  when  fluid  pareels  exehange  energy  vertieally  by  eollision.  The  transport  terms 
only  move  the  total  energy  around  vertieally  and  they  never  ehange  it.  Consequently  if 
the  sign  of  these  terms  is  negative  at  a  given  height,  it  must  be  positive  at  another  height. 


The  fourth,  or  e  terms  are  the  viseous  dissipation  terms  and  are  effeetive  at  the 
smallest  seale  (Kolmogorov  miero  seale)  of  the  turbulenee  where  moleeular  viseosity 
overwhelms  inertial  eomponents  of  the  turbulenee. 

The  fifth  terms  in  eaeh  eomponent  are  pressure-strain  eorrelation  terms.  These 
terms  explain  the  energy  transfer  among  the  three  eomponents.  Pressure  ean  also  transfer 
the  kinetie  energy  from  one  direetion  to  the  lateral  direetion.  Sinee  the  TKE  is 
eonservative,  the  summation  of  these  terms  in  eaeh  eomponent  will  be  zero. 

The  last  terms  represents  the  Coriolis  redistribution  of  a  fluid  pareel.  The  motion 
of  a  fluid  pareel  ehanges  its  direetion  under  the  Coriolis  foree  to  the  right  in  northern 
hemisphere  and  to  the  left  in  southern  hemisphere,  but  the  total  kinetie  energy  remains 
the  same.  The  summation  of  these  terms  in  eaeh  eomponent  is  zero. 


The  total  TKE  equation  ean  be  obtained  by  adding  up  Equations  2.29,  2.30  and 
(2.31).  The  total  TKE  is. 


_a 

dt 


,2  ,2  ,2 

u  -\-v  -\-w 
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du  -r^  dv 

—u  w - vw  — 

dz  dz 
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z' 

,2  ,2  ,2  , 
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(2.32) 


The  dissipation  term  has  nine  eomponents  in  it,  whieh  are  the  eombination  of  eaeh 
TKE  eomponent’ s  dissipation  terms. 


£  =  £  ,2  +£  .2  +£  ,2 


=  u 


ydx  j 


ydXj 


dx 


^  du^ 
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du  ) 

,  f ^  J 

^  dv^ 
[^yj 

+  I 

'  dv'' 
^dz) 

'  dw'\ 

^  dw 

V 

(2.33) 
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At  the  finest  scales,  Kolmogorov  micro  scale,  turbulence  is  nearly  isotropic  and  there  is 
no  direction  preference  at  these  scales,  so  each  component  of  s  is  1/3  of  the  total 
dissipation.  Furthermore,  the  total  dissipation  is  approximately  fifteen  times  any  one  of 
the  diagonal  terms  in  (2.33)  due  to  the  isotropy.  Hence  s  can  be  presented  as. 


(2.34) 


The  energy  spectrum  of  the  total  TKE  shows  us  how  strong  the  turbulent  velocity 
is  at  different  time  or  length  scales.  The  log-log  plot  of  the  wave  number  and  the 
correspondent  energy  at  each  wave  number  for  a  steadily  forced  flow  is  shown  in  Figure 
(2.2). 


Figure  2.2  Energy  Density  Spectrum  of  TKE;  k  is  wave  number  (1/cm),  E  (k)  is  the 
Energy  as  a  function  of  k 

Figure  2.2  shows  the  energy  density  spectrum  of  TKE  for  a  steadily  forced 
turbulent  fluid  with  forcing  scales  much  larger  than  dissipation  scales.  If  the  Reynolds 
number  is  much  larger  than  one,  which  means  small  molecular  viscosity,  an  inertial-sub 
range  occurs.  The  wider  the  (in  wave  number  space)  inertial  sub  range,  the  larger  the 
TKE.  If  the  slope  of  the  inertial  sub  range  is  at  least  ,  there  is  a  turbulent  flow. 
Energy  enters  the  flow  at  low  frequencies,  at  macro  scale,  and  then  this  energy  is 
transported  from  larger  eddies  to  the  smaller  eddies  along  the  energy  density  spectrum. 
A  is  the  integral  or  macro  scale  of  turbulent  flow.  This  is  the  scale  where  the  largest 
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eddies  oeeur.  rj  is  the  Kolmogorov  or  miero  seale  of  the  turbulent  flow  and  is  the  seale  of 
the  smallest  eddies  where  moleeular  viseosity  beeomes  important,  and  finally  dissipation 
of  the  turbulent  flow  oeeurs. 

The  basie  equations  related  to  turbulent  flow  and  their  physieal  meanings  were 
explained  sinee  the  Langmuir  eireulations,  the  main  foeus  of  this  study,  are  eoherent 
motions  that  eo-exist  with  turbulent  eddies  in  oeean  mixed  layers  in  the  presenee  of  wind 
foreing  and  waves.  These  eireulations  play  an  important  role  for  the  mixing  in  the  upper 
oeean.  A  review  of  these  eireulations  will  be  given  in  the  next  seetion. 


B,  LANGMUIR  CIRCULATIONS 

Langmuir  eireulations  are  parallel  series  of  the  eounter  rotating  helieal  vortiees  in 
the  surfaee  layer  of  the  oeeans  and  the  lakes  having  their  axes  aligned  within  a  few 
degrees  along  the  wind  direetion  in  the  presenee  of  surfaee  gravity  waves.  Beeause  of 
their  long-lived,  signifieantly  large-seale  eireulation  within  the  oeean  mixed  layer,  it  is 
hypothesized  that  they  eontribute  to  the  heat  and  momentum  transfer  in  the  mixing  layer 
in  oeeans  and  lakes. 


Figure  2.3  The  colleetion  of  the  bubbles  into  long  and  narrow  streaks  on  the  ocean 
surface 

Langmuir  (1938)  showed  that  the  collection  of  the  seaweed,  marine  organisms, 
and  foam  from  breaking  waves  or  oil  leakages  into  streaks,  like  those  in  Figure  2.3,  is  the 
visible  evidence  of  counter  rotating  helical  vortices  near  the  surface.  These  features  are 
known  as  Langmuir  circulations  after  Dr.  Irving  Langmuir  who  was  the  first  scientist  to 


13 


study  them  and  defined  their  properties  without  explaining  the  formation  meehanisms. 
The  most  significant  features  of  Langmuir  circulations  are  lateral  coherence  and  long  life, 
which  makes  them  different  from  eddies  typical  of  turbulent  flow. 

The  sketch  in  Figure  2.4  shows  the  basic  structure  of  the  Langmuir  circulations  as 
described  by  Langmuir  (1938)  and  following  investigators. 


Figure  2.4  A  schematic  of  forcing  factors  and  structures  of  Langmuir  circulations 
(From  Smith  2001) 

In  Figure  2.4,  L  is  the  horizontal  length  in  the  crosswind  direction  between  two 
convergence  zones  or  two  divergence  zones,  also  known  as  row  spacing.  D  is  the  vertical 
length  scale  or  penetration  depth,  and  is  often  limited  by  the  pycnocline  depth  (Langmuir 
1938).  The  common  coordinate  system  used  by  most  of  the  investigators  is  the  wind- 
based  coordinate  system,  which  represents  x  in  along  wind  direction,  z  in  upward  and  y  is 
cross  wind  direction.  Velocity  components  u,  v  and  w  are  in  the  directions  of  x,  y  and  z 
respectively.  Wind  speed  {U^^,  generally  defined  10  m  above  the  water  surface,  is  also 

an  important  parameter  for  the  Langmuir  circulations  since  wind  is  believed  to  be  a  main 
factor  forming  these  circulations.  Observations  show  that  there  is  downwelling  under  the 
convergence  zones  and  upwelling  under  the  divergence  zones.  The  u  velocity 
component  is  maximum  over  the  convergence  zones  and  minimum  over  the  divergence 
zones  as  shown  in  Figure  2.4. 
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1,  Formation  Theories  of  Langmuir  Circulations 

There  have  been  many  theories  proposed  for  the  formation  of  Langmuir 
eireulations  where  the  observed  features  of  the  cireulations  are  attempted  to  be  deseribed 
by  physieal  models.  However  there  is  still  not  a  theory,  whieh  ean  explain  all  the 
observed  features  exaetly.  The  most  commonly  accepted  theories  are  related  to  the  wave- 
current  interactions.  A  review  of  Langmuir  circulation  theories  follows. 

Thermal  convection  has  been  hypothesized  to  be  main  formation  mechanism  by 
several  investigators  following  observations  of  collections  of  convective  plumes  into 
rows  by  the  wind  shear  in  the  atmosphere  (Pollard  1977).  However,  Langmuir  (1938) 
advocated  that  the  wind  is  the  primary  mechanism  for  the  circulations,  supported  by 
observed  circulations  under  stable  or  neutral  conditions.  Stommel  (1951),  Faller  (1964), 
Csanady  (1965),  Ichye  (1967)  and  Myer  (1969)  have  all  observed  the  Langmuir 
circulations  in  conditions  of  stable  stratifications  (Pollard  1977).  Li  and  Garrett  (1995) 
found  that  wave  forcing  dominates  over  the  thermal  convection  for  typical  wind  speed 
and  surface  cooling  in  driving  Langmuir  circulations.  Since  Langmuir  circulations  are 
observed  under  stable  conditions,  thermal  convection  cannot  be  considered  as  a  primary 
formation  mechanism,  but  may  have  secondary  effects  contributing  to  the  primary 
mechanism.  For  example,  under  the  surface  cooling  conditions,  circulations  may  occur 
with  less  wind  speed  than  the  stable  conditions,  while  stratification  induced  by  strong 
surface  heating  may  suppress  the  circulations  generated  by  the  wave  forcing  (Li  and 
Garrett  1995). 

Coupling  of  the  ocean  with  the  atmospheric  rolls  has  also  been  considered  as  a 
forcing  mechanism  for  the  circulations.  However  the  air  moving  with  these  rolls  advects 
too  quickly  across  the  surface  to  couple  with  the  water  and  has  negligible  effect  on  the 
water.  Pollard  (1977)  cited  that  the  observations  of  sea  smoke  above  lakes  made  by 
Langmuir  (1938)  and  Stommel  (1951)  support  this  conclusion.  Leibovich  (1983)  also 
stated  that  the  position  of  these  rolls  do  not  align  with  the  windrow  streaks  above  the 
Langmuir  circulations. 

Damping  of  the  capillary  waves  in  slicks  occurs  in  the  presence  of  surface  films. 
Welander  (1963)  and  Kraus  (1967)  proposed  that  the  effects  of  these  surface  films  might 
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be  the  main  process,  which  create  the  Langmuir  circulations.  Welander  (1963) 
hypothesized  that  smoothing  the  surface  by  damping  the  capillary  waves  in  slicks  would 
cause  a  higher  wind  speed,  accelerating  the  water,  then  producing  a  secondary  inflow  into 
the  slicks  region.  Leibovich  (1983)  discounted  this  suggestion  since  increasing  wind 
speed  is  the  result  of  the  decreasing  surface  drag,  and  therefore  the  water  speed  in  slicks 
would  be  smaller  than  the  surrounding  water,  which  is  not  observed.  Kraus  (1967) 
suggested  another  surface  film  related  mechanism  where  energy  is  transferred  from 
capillary  waves  into  the  slicks  while  they  are  being  damped  as  they  approach  the  slicks 
accelerating  the  slicks  in  the  wind  direction.  Pollard  (1977)  argues  that  it  is  doubtful 
whether  the  capillary  waves  can  supply  energy  fast  enough  to  be  sufficient  for  the 
observed  rapid  generation  of  Langmuir  circulation.  The  energy  budget  given  by  Myer 
(1971)  supports  Pollard’s  argument  since  the  total  energy  transfer  from  the  capillary 
waves  to  the  slicks  is  about  10  %  of  the  rate  needed  to  generate  the  circulations. 
Furthermore  Leibovich  (1983)  stated  that  the  formation  of  circulations  in  laboratory 
experiments  exists  in  the  absence  of  surface  films. 

Ekman  layer  instabilities  may  result  as  convective  rolls  with  a  small  angle  to  the 
wind,  and  Faller  (1964)  suggested  that  these  instabilities  of  Ekman  spiral  could  be  the 
driving  mechanisms  for  Eangmuir  circulations.  Eeibovich  (1983)  stated  that  the  growth 
rates  of  such  instabilities  are  too  small  to  account  for  Eangmuir  circulations.  Gammelsrod 
(1975)  enhanced  the  Ealler  (1964)  idea  by  starting  with  a  steady  along  wind  shear  flow, 
which  is  the  result  of  the  balance  between  Coriolis  force  and  friction.  The  vertical 
advection  of  the  mean  shear  by  vertical  perturbation  velocity  will  accelerate  the  along 
wind  perturbation  flow  and  then  the  resultant  along  wind  velocities  will  be  turned  to  the 
cross  wind  direction  by  the  Coriolis  force.  As  a  result,  the  Eangmuir  like  cellular 
structures  will  occur  due  to  the  new  cross  wind  shear  flow.  Gammelsrod’s  theory  predicts 
the  cell  orientation  at  an  angle  to  the  vertical  and  along  wind  velocities  at  the  surface  to 
be  zero,  which  are  the  main  differences  from  the  observed  vertical  Eangmuir  cells. 
Eurthermore,  the  Coriolis  force  acting  on  the  along  wind  perturbation  velocities  is  not 
balanced  by  friction  and  remain  long  enough  to  create  the  crosswind  velocities.  This  is 
not  realistic,  because  the  friction  is  created  by  small-scale  turbulence  via  Reynolds 
stresses  and  this  feature  is  displaced  vertically  by  itself  and  the  friction  term  will  continue 
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to  balance  the  Coriolis  force  on  the  perturbation  horizontal  veloeities  preventing  the 
growth  of  erosswind  veloeities  (Pollard  1977). 

Theories  based  on  wind  foreing  meehanisms  in  the  presenee  of  surfaee  gravity 
waves  are  the  most  eonsistent  with  the  observations.  All  the  experimental  studies 
ineluding  Langmuir’s  own  indieate  that  the  wind  is  the  primary  driving  mechanism  for 
the  Langmuir  eireulations.  However,  the  question  of  how  the  wind  sets  up  these 
eireulations  is  still  not  elear.  The  wind-based  theories  ean  be  generalized  as  waves  and 
wind  drift  eurrent  interaetion  theories.  Craik  and  Leibovieh  (1976)  and  Garrett  (1976)  are 
two  leading  studies  that  propose  the  most  aeeepted  wave  and  eurrent  interaetion  theories. 
Both  theories  assume  a  nearly  irrotational  wave  field  and  a  mueh  weaker  rotational 
eurrent  in  upper  layer,  with  the  nonlinear  interaetions  between  the  waves  and  eurrents 
forming  Langmuir  eireulations.  The  meehanism  of  how  the  interaetions  ereate  the 
eireulations  is  expressed  differently  in  these  two  theories. 


Figure  2.5  Refraetion  of  a  wave  ray  by  a  eurrent  U  (y).  The  arrows  indieate  the  lateral 
forees  exerted  by  the  wave  (From  Garrett  1976) 

Garrett  (1976)  related  Langmuir  eireulations  to  the  refraetion  of  the  surface  waves 
by  the  mean  eurrent  anomaly,  whieh  is  in  the  direetion  of  the  winds.  Figure  2.5  illustrates 
the  refraetion  of  the  waves  by  the  mean  eurrent  jet  U(y).  The  first  meehanism  of  this 
refraetion  is  the  lateral  forees  that  tend  to  produee  an  inflow  towards  the  maximum 
downwind  eurrent.  The  wave  speed  perpendieular  the  mean  eurrent  will  deerease  in  the 
current  zone  and  then  inerease  out  of  the  eurrent  zone,  so  there  must  be  a  momentum 
transfer  between  the  waves  and  eurrents  to  aeeomplish  this  aetion,  and  the  lateral  forees 
are  the  result  of  this  momentum  transfer.  In  Figure  2.6,  the  arrows  indieate  the  forees 
exerted  by  waves  propagating  in  positive  x-y  direetions.  Garrett  showed  that  the  foree  is 
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predominantly  towards  the  eurrent  maximum  for  waves  at  a  small  angle  to  the  x-axis. 
The  inflow  produeed  by  dominant  lateral  forees  ereates  along  wind  vortieity,  whieh  is  the 
signature  of  Langmuir  eireulations. 


y 


Figure  2.6  Lateral  forees  produeed  by  refraetion  of  the  waves  (From  Garrett  1976) 

Enhaneed  wave  breaking  is  the  seeond  eonsequenee  of  the  wave  refraetion  effeet. 
(Garrett  1976).  The  wave  amplitude  will  inerease  as  the  wave  speed  is  deereasing  in  the 
high  eurrent  zone  and  this  will  inerease  the  wave  breaking  in  that  region  resulting  the 
momentum  transfer  from  the  waves  to  the  eurrents  whieh  will  reinforee  U(y).  The 
momentum  transfer  in  the  y-direetion  is  eaneeled  sinee  the  wave  speetrum  is  assumed  to 
be  symmetrie  about  the  downwind  direetion  (Garrett  1976).  Enhaneed  U(y)  ereates  a 
feedbaek  for  the  first  meehanism,  thus  the  eombination  of  two  meehanisms  produees  an 
instability  meehanism,  whieh  forees  the  Eangmuir  eireulations. 

Eigure  2.7  shows  that  the  vortex  line  of  z-vortieity  is  tilted  towards  the  downwind 
direetion  by  the  Stokes  drift  of  the  waves.  Aeeording  to  the  Garret’s  theory,  the  z- 
vortieity  with  a  vertieal  vortex  line  is  ereated  by  the  wave  dissipation  in  the  y  direetion 
and  the  vortex  line  is  bent  in  the  x  direetion  by  the  x-eomponent  of  the  Stokes  drift. 
Einally  the  eellular  motion  will  be  about  this  new  x-oriented  vortex  line  ereating  the  x- 
vortieity  as  a  signal  of  the  Eangmuir  eireulations. 
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Figure  2.7  The  z-  vortieity  beeomes  x-  votieity  due  to  tilting  of  the  z-  oriented  vortex 
line  towards  the  direetion  of  maximum  U(y)  by  the  Stokes  drift  (From 
Garrett  1976) 

Garrett  (1976)  also  stated  the  basie  limitations  of  his  theory.  The  internal 
refleetion  of  the  waves  by  the  eurrent  is  negleeted.  For  the  small  wave  angles  to  the 
eurrent  zone  used  in  the  theory,  it  has  been  found  that  the  weak  eurrents  ean  internally 
refleet  the  waves  whieh  ereates  the  inflow  towards  the  line  of  refleetion  instead  of 
maximum  U(y)  region.  Although,  the  eonvergenee  zone  exists  in  a  different  plaee  than 
Garrett  (1976)  explained  previously,  Garrett  (1976)  asserts  that  lateral  forees  are  still 
exerted  by  the  refleeted  waves  and  the  feedbaek  meehanism  for  the  instability  remains 
the  same.  The  assumption  of  slow  variation  of  U(y)  was  stated  as  another  limitation.  The 
wave  modifieation  by  the  eurrents  is  mathematieally  and  physieally  explained  by 
assuming  mean  steady  flow  while  the  mean  flow  is  being  enhaneed  in  time,  and  there 
seems  to  be  a  eontradietion  between  the  assumption  and  the  applieation.  Although  Garrett 
(1976)  explained  this  ineonsisteney  as  a  limitation,  he  mentioned  that  the  developing  time 
seale  of  the  mean  flow  is  mueh  slower  than  the  set  up  time  seale  of  the  Langmuir 
eireulations  by  the  wave  paeket.  Another  possible  limitation  given  by  Garrett  (1976)  is 
the  eonstant  surfaee  eurrents  over  the  depth  of  waves  and  he  emphasized  that  even  though 
the  eurrent  is  not  uniform  over  this  depth;  the  interaetion  stresses  exerted  by  the  waves 
will  still  eause  the  vertieal  shear. 


Leibovieh  (1983)  stated  that  the  modifieation  of  the  wave  field  by  the  horizontal 
eurrent  and  the  preferential  wave  breaking  are  doubtful  theoretieally  and  experimentally. 

Leibovieh  (1983)  elaimed  that  WKBJ  theory  used  by  Garrett  to  explain  the  amplifieation 

19 


of  the  waves  and  may  lead  to  wrong  eonclusions.  Leibovieh  (1983)  eited  that  Smith 
(1980)  was  unable  to  find  amplifieation  of  the  waves  predieted  by  Garrett’s  (1976)  model 
without  WKBJ  theory.  Leibovieh  (1983)  also  eited  the  experimental  studies  done  by 
Thorpe  and  Hall  (1980)  and  Kenney  (1977),  and  these  experimental  studies  found  the 
wave  breaking  between  streaks  and  in  them  with  equal  frequeney  instead  of  dominantly 
on  the  line  of  eonvergenee. 

Pollard  (1977)  suggested  that  Garrett’s  theory  is  plausible  sinee  it  ean 
qualitatively  explain  all  observed  features  but  he  eritieized  the  mathematies  of  this  theory 
in  terms  of  the  same  aspeets  as  Leibovieh  (1983) 

The  other  important  theory  related  to  the  wave  and  eurrent  interaetions  is  the 
Craik  and  Leibovieh  (CL)  (1976)  theory.  Craik  and  Leibovieh  (1976)  defined  the 
eommon  qualitative  observed  features  of  the  eireulations  that  a  theory  should  be  able  to 
explain: 

CL  (1)  A  parallel  system  of  vortiees  aligned  with  the  wind  must  be  predieted. 

CL  (2)  A  means  must  be  given  by  whieh  these  vortiees  are  driven  by  the  wind. 

CL  (3)  The  resulting  eells  must  have  the  possibility  of  an  asymmetrie  strueture 
with  downwelling  speeds  larger  than  the  upwelling  speeds. 

CL  (4)  Downwelling  must  be  below  where  the  wind  direeted  surfaee  eurrent  is 
greatest 

CL  (5)  The  Langmuir  eireulations  must  have  maximum  downwelling  speeds 
comparable  to  the  mean  wind  directed  surface  currents. 

There  are  two  distinct  versions  of  CL  theories.  The  first  one,  CL-1,  was  developed 
by  Craik  and  Leibovieh  (1976)  and  the  second  one,  CL-2,  was  first  proposed  by  Craik 
(1977)  and  then  explored  by  Leibovieh  (1977a- 1977b).  Both  versions  depend  upon  the 
discussion  of  the  “Vortex  Force”,  which  is  the  cross  product  of  Stokes  drift  and  the  mean 
vorticity, 

f -u^xco  (2.35) 


20 


where  (o  is  the  mean  vortieity  and  is  the  Stokes  drift  of  surfaee  waves.  Aeeording  to 

the  small  amplitude  sinusoidal  wave  theory,  water  partieles  move  in  elosed  eireular 
orbits,  and  there  is  no  net  motion  of  the  water  mass.  However  in  real  water  waves  of 
finite  amplitude,  the  orbits  are  not  elosed,  therefore  there  is  a  small  amount  of  net  motion 
of  water  mass  in  the  wave  traveling  direetion.  This  is  ealled  Stokes  drift  and  Philips 

(1977)  defined  the  Stokes  drift  as  where  braekets  represent  the 

temporal  mean,  w^is  the  orbital  veloeity  of  the  water  mass.  Basieally,  CL  theories  are 
based  on  the  restoring  of  the  mean  vortieity  by  surfaee  wave  Stokes  drift  via  vortex  foree. 


Craik  and  Leibovieh  (1976)  introdueed  the  CL-1  theory  by  evolving  the 
previously  developed  equations  by  Craik  (1970)  and  Leibovieh  and  Ulrieh  (1972).  The 
initial  assumptions  are  the  eoverage  of  the  surfaee  by  short  erested  wind  generated 
irrotational  gravity  waves,  whieh  have  a  diserete  direetional  frequeney  speetrum  with  two 
peaks  at  an  equal  angle  6  in  both  sides  of  the  wind  direetion  (monoehromatie  bimodal 
wave),  infinite  depth,  eonstant  eddy  viseosity  to  represents  the  turbulenee  effeets,  and  a 
wave  spatial  strueture  that  remains  eoherent  for  Stokes  drift  averaging.  Craik  and 
Leibovieh  (1976)  explained  the  Stokes  drift  is  periodie  in  the  y  direetion  (eross-wind 
direetion)  for  the  monoehromatie  bimodal  wave  assumption. 
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Figure  2.8  CL  -  1  instability  mechanism  creating  Langmuir  circulations.  C  is  the 
vortex  force,  a  is  y-vorticity,  and  Us  is  the  Stokes  drift  in  the  x-direction. 
The  coordinate  system  is  the  wind-based  coordinate  system  with  the 
positive  x-axis  towards  the  wind  direction.  (From  Leibovich  1983) 

Leibovich  (1983)  illustrated  the  CL-1  mechanism  (Figure  2.8).  The  coordinate 

system  in  Figure  2.8  is  the  wind-based  coordinate  system.  When  the  wind  blows  in  a 

fixed  direction,  the  initial  wave  field  will  have  a  unidirectional  Stokes  drift  in  the  same 

direction  as  the  wind.  This  initial  Stokes  drift  is  =  UJ ,  where  is  the  total  Stokes 

drift  magnitude  (Leibovich  1983).  The  momentum  transfer  from  the  wind  to  the  water 
dictates  the  horizontal  velocity  in  the  same  direction  with  the  wind.  Horizontally 
averaged  velocity  is  defined  as  u  =  f/(z,t)z  with  positive  y-vorticity  a)  =  {dU ldz)j  by 

neglecting  the  Coriolis  force  (Leibovich  1983).  The  cross  product  between  the  Stokes 
drift  and  the  y-vorticity  using  Equation  2.35  will  give  an  upward  vortex  force. 

l=kU^{dU/dz)  (2.36) 

Leibovich  (1983)  explained  that  this  vortex  force  is  analogous  to  the  buoyancy 
force,  and  if  the  horizontal  velocity  is  only  a  function  of  depth,  the  vortex  force  will  be 
balanced  by  the  hydrostatic  pressure  gradient,  which  produces  unidirectional  current. 
After  the  characteristic  of  initial  wave  field  turns  into  the  assumed  wave  field,  the  Stokes 
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drift  will  be  periodic  in  the  y-direction.  The  total  Stokes  drift,  periodic  in  y-direction  is 
given  by  Craik  and  Leibovich  (1976)  in  Equation  2.37. 

=  2ke^^  [l  +  cos  2/y ]  (2.37) 

Where  k  =  cos0  and  /  =  sin^  are  the  wave  numbers  in  x-  and  y-directions  respectively. 
0  is  the  angle  between  the  wind  direction  and  the  crests  of  the  bimodal  short-crested 
surface  gravity  waves  on  each  side  of  the  wind.  Since  the  Stokes  drift  will  vary  in  the  y- 
direction,  the  upward  vortex  force  will  not  be  balanced  by  the  hydrostatic  pressure 
gradient  anymore  and  will  also  change  in  the  y-direction  as  a  result  of  rectified  Stokes 
drift.  The  variations  of  the  upward  vortex  force  in  the  y-direction  will  produce 
overturning  about  the  x-axis.  In  Figure  2.8,  the  Stokes  drift  at  the  wave  crests  is  higher 

than  the  Stokes  drift  between  the  wave  crests,  and  the  upward  vortex  force  is  maximum 
and  minimum  at  the  same  places  where  the  Stokes  drift  is  maximum  and  minimum 
respectively.  The  variation  of  the  upward  vortex  force  in  the  y-direction  produces  two 
circular  counter-rotating  rolls  converging  at  the  surface  between  the  wave  crests.  These 
two  counter-rotating  circular  motions  are  the  Langmuir  circulations  as  described  by  the 
CL-1  theory. 

Craik  and  Leibovich  (1976)  reported  that  their  model  results  qualitatively 
matched  with  the  observed  circulation  characteristics.  Nevertheless,  for  realistic  wind 
speeds  and  the  wave  slopes,  Craik  and  Leibovich  (1976)  were  forced  to  admit  that  their 
model  failed  to  satisfy  requirement  CL  (3).  Pollard  (1977)  noted  that  the  CL-1  model 
forecasts  upwelling  regions  where  the  Stokes  drift  is  largest,  which  is  inconsistent  with 
observations.  Pollard  (1977)  also  emphasized  that  such  a  Stokes  drift  pattern  as  explained 
in  the  CL-1  theory  requires  that  variances  in  the  y-direction  should  remain  fixed  in  space 
for  a  long  enough  time  to  produce  Langmuir  circulations,  and  crossed  wave  pairs  should 
be  phase  locked  for  several  wave  periods  to  achieve  this.  This  type  of  phase  locking  is 
not  expected  in  a  wind-generated  sea,  although  it  can  be  generated  in  laboratory 
experiments  (Leibovich  1983). 

It  is  generally  accepted  that  the  assumptions  of  time  independent  currents,  crossed 
waves  and  invariance  of  rectified  motion  in  the  x-direction  are  the  weakest  points  of  the 
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CL-1  theory.  Leibovich  (1977a)  modified  the  theory  to  inelude  time  evolution  of  the 
wind-direeted  eurrents  and  eireulations  under  the  assumption  that  the  wave  field  was 
symmetrie  with  respeet  to  the  x-axis  and  invariant  in  the  x-direetion.  This  assumption 
made  the  motion  independent  of  x,  and  the  formulation  of  the  phenomena  was  explained 
as  an  initial  value  problem  where  the  eurrents  and  eireulations  were  initially  zero.  The 
dimensionless  Langmuir  number,  La,  was  derived  by  Leibovieh  (1977a)  to  represent  the 
initial  value  problem. 

La-[olk^ ! (ja^ul^  (2.38) 

In  this  expression,  is  the  eddy  viseosity,  cr  is  the  surfaee  waves  frequeney,  m,  is  the 

frietion  veloeity  for  eonstant  wind  stress  when  t  >  0 ,  k  is  the  wave  number,  and  a  is  the 
wave  amplitude.  The  Langmuir  number,  whieh  is  the  simply  the  inverse  of  the  Reynolds 
number,  physieally  expresses  the  balanee  between  the  rate  of  diffusion  and  the  rate  of 
produetion  of  stream-wise  vortieity  (Leibovieh  1983). 

Leibovieh  (1977a)  also  assumed  infinite  depth,  eonstant  eddy  viseosity,  eonstant 
density,  and  that  the  Stokes  drift  was  produced  by  a  pair  of  uniform,  narrow-banded  wave 
trains  propagating  at  equal  and  opposite  angles  with  respect  to  the  wind  direction. 
Leibovich  (1977a)  used  the  Langmuir  number  to  define  the  scales  of  Langmuir 
circulations,  as  shown  in  Table  1. 


Length 

Time(7;) 

Wind-directed 

component  of 

rectified  current 

Crosswind  velocity 

component  of 

rectified  currents 

/ <j^'^  jaku-t 

ut  jvjk 

Table  1.  Scales  of  Langmuir  circulations  defined  by  Leibovich  (1976a)  (From 
Leibovich  1983) 


The  improved  version  of  the  CL-1  mechanism  by  Leibovich  (1977a)  satisfies  the 
listed  features  CL  (1)  through  CL  (5)  and  matches  observed  results  better  than  its 
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predecessor.  In  the  revised  model,  the  maximum  surface  sweeping  speeds  are  close  to  the 
maximum  downwelling  speeds,  and  the  maximum  downwelling  speeds  range  from 
slightly  less  than  twice  to  more  than  seven  times  the  maximum  upwelling  speeds.  The 
maximum  along-wind  current  anomaly  is  found  over  the  downwelling  zones  with  a  speed 
approximately  55-75  %  greater  than  the  speeds  midway  between  zones.  (Leibovich 
1983). 

Craik  (1977)  first  introduced  the  second  version  of  the  CL  theory  (CL-2)  for  a 
constant  density  system.  Leibovich  (1977b)  adapted  the  theory  to  include  density 
stratifications  under  the  Boussinesq  approximation.  The  time-dependent  development  of 
the  waves  by  the  wind  is  neglected  in  CL-2  theory,  and  both  momentum  and  thermal 
eddy  diffusivities  are  assumed  to  be  constant.  It  is  necessary  to  have  a  non-coherent 
surface  wave  structure  (Leibovich  1983).  The  vortex  force  is  again  the  reason  for  the 
circulations,  but  the  CL-2  mechanism  producing  the  vortex  force  is  different  from  that 
used  by  the  CL-1  theory.  Leibovich  (1983)  illustrates  both  the  dynamics  and  kinematics 
of  the  CL-2  mechanism,  as  shown  in  Figure  2.9. 


Figure  2.9  CL-2  instability  mechanism  of  Langmuir  circulation  generation.  C  is  the 
vortex  force,  (o  is  the  z-vorticity,  and  Us  is  the  Stokes  drift  in  the  x- 
direction  The  coordinate  system  is  the  wind-based  coordinate  system  with 
the  positive  x-axis  towards  the  wind  direction.  (From  Leibovich  1983) 
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Stokes  drift  is  defined  as  u^=U^[z)i  if  the  waves  lack  a  coherent  spatial 
structure.  Decreasing  and  U  with  depth  will  lead  to  decreasing  vortex  force  with 
depth,  and  combined  results  of  such  a  distribution  of  and  U  are  similar  to  those 

obtained  with  a  statically  unstable  density  distribution  (Leibovich  1983).  If  dissipation  is 
weak  enough,  the  current  is  unstable.  Horizontal  uniform  current  U (z)  is  unidirectional 
in  the  wind  direction,  because  the  vortex  force  in  Equation  2.35  is  balanced  by  a  pressure 
gradient  force  in  the  absence  of  the  disturbances  in  the  y-direction,  as  assumed  in  CL-1 
theory.  The  vertical  vorticity  co^  is  produced  by  the  irregularities,  u{y,z,t) ,  in  the  largely 

horizontally  uniform  current  U{z).  Equation  2.5  becomes  co^  =-{duldy)k  (  where  k  is 
the  unit  vertical  vector).  The  vortex  force  per  unit  mass  is  then  defined  as  in  Equation 
2.39. 

f.=-U,oyJ  2.39 

If  the  z-vorticity  co^  is  positive  (upward),  vortex  force  will  be  in  the  negative  y- 
direction,  and  if  z-vorticity  co^  is  negative  (downward),  vortex  force  will  be  in  the  positive 

y-direction.  In  both  cases  the  water  mass  accelerates  towards  the  maximum  current 
anomaly,  so  convergence  zones  will  be  created  in  this  region.  Because  of  the 
conservation  of  mass,  water  at  the  convergence  zones  will  sink,  and  the  sinking  water 
will  be  replaced  by  the  horizontally  converging  flow.  This  increases  the  magnitude  of  the 
convergent  jet  due  to  the  vortex  force.  In  summary,  the  rotation  and  stretching  of  the 
vertical  vorticity  by  Stokes  drift  results  in  convergence  and  amplification  of  the  surface 
current  anomaly.  This  instability  process  is  the  basis  of  CE-2  theory  and  creates  the 
Eangmuir  cells. 

Eeibovich  (1977b)  estimated  the  instability,  M(z),  of  the  system  condition  for  an 
inviscid  stratified  fluid  in  infinite  depth. 

M{z)  =  U[{z)U'{z)-N\z)  (2.40) 

The  instability  occurs  when  M{z)  is  positive  and  does  not  occur  when  M(z)  is 
negative.  N^{z)  is  the  Brunt- Vaisala,  or  buoyancy,  frequency  of  the  initial  state  as 
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described  in  Equation  2.2.  However,  it  can  be  also  formulated  in  terms  of  thermal 
expansion  as  /3gT  (z) ,  where  P  is  the  thermal  expansion  coefficient  and  T  (z)  is  the 
vertical  gradient  of  temperature  distribution.  U\  is  the  vertical  gradient  of  the  Stokes  drift 
and  U  is  the  vertical  gradient  of  the  shear  currents.  M(z)  is  maximum  at  the  surface 
(z  =  0).  Leibovich  (1977b)  defined  the  maximum  growth  rate  of  the  circulations  for 
unstable  system  conditions  with  infinitesimal  wavelengths  as 

(2-41) 

When  the  instability  of  the  system  occurs  with  stable  density  stratification,  there 
will  be  a  depth  below  which  the  disturbances  will  not  penetrate.  Leibovich  (1983) 
expressed  the  instability  condition  in  terms  of  the  Richardson  number,  Ri. 

Ri  =  mm[N^lu'{z)U[{z)\<\  (2.42) 

Leibovich  (1983)  cited  the  studies  Leibovich  and  Paolucci  (1980a,  b,  1981)  and 
Leibovich  and  LeLe  (1982),  which  have  explored  the  CL-2  mechanism  with  nonzero 
eddy  diffusivities.  Constant  positive  and  a  time-dependent  unidirectional  current 
solution  of  the  equations  of  motion  were  considered  in  Leibovich  (1977b)  for  infinitely 
deep  water  and  constant  surface  stress  at  t  =  0.  Leibovich  and  Paolucci  (1980a)  found 
that  the  small-scale,  very  weak  disturbances  grew  in  strength  while  they  cascaded  to 
larger  scales  unless  the  domain  size  prevented  them  from  growing.  Leibovich  and 
Paolucci  (1981)  calculated  the  linear  stability  limits  for  two-dimensional  rolls  in  terms  of 
Richardson  number,  inverse  Langmuir  Number  and  length  scale  (Table  2). 


Ri 

La: 

K 

0.00 

1.52 

0.32 

0.10 

1.58 

0.31 

0.25 

1.66 

0.30 

Table  2.  Critical  inverse  Langmuir  number  (  La^ ' )  and  critical  wave  number  ( )  of  linear 

theory  at  various  Richardson  numbers  {Ri)  defined  by  Leibovich  and  Paolucci 
(1981)  (Prom  Leibovich  1983) 
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Leibovich  and  Paolucci  (1980a)  expressed  the  Richardson  number  in  Table  2 
with  Equation  2.43  for  a  temperature-stratified  ocean 


Ri  =  /3gT  (^au^ky  <j I Oj. 


(2.43) 


where  k  is  the  dominant  wave  number.  The  convective  motions  mix  the  upper  layers  and 
form  a  layer  of  strong  temperature  gradient.  After  correcting  several  errors  in  the  work  of 
Leibovich  and  Paolucci  (1980a),  Leibovich  and  LeLe  (1982)  found  smaller  penetration 
depth  of  Langmuir  circulations.  They  also  found  that  preexisting  thermoclines  act  like  a 
strong  obstacle,  which  prevent  convective  motions  from  deeper  penetrations  (Leibovich 
1983). 

Leibovich  (1983)  recommended  that  both  CL  theories  be  tested  for  conditions  in 
which  waves  and  currents  exist  in  the  absence  of  the  wind,  such  that  is  not  parallel  to 
the  horizontally  averaged  current  U .  The  CL  theories  suggest  that  an  instability  will 
occur  whenever  u^U  >  0 ,  whether  this  results  in  Langmuir  circulations  or  not.  Leibovich 
(1983)  also  noted  that  the  CL  theory  of  Langmuir  circulations  in  two  dimensions  is 
analogous  to  the  double  diffusive  convection  problems  mathematically,  but  Leibovich 
and  Tandon  (1992)  stated  that  this  is  not  true  for  CL  theory  in  three  dimensions. 

Craik  (1982a)  proved  that  water  waves  induce  the  weak  viscous  Eulerian  current 
by  the  instability  of  CL-2  theory  in  the  absence  of  wind,  and  noted  that  this  may  be  the 
same  case  for  the  currents  induced  by  internal  waves  near  a  density  interface.  Leibovich 
(1983)  cited  this  study  as  one  to  be  explored  in  terms  of  the  CL  theories. 

After  reviewing  previous  analysis,  numerical  models,  laboratory  simulations,  and 
filed  work.  Smith  (2001)  concluded  that  the  CL-2  mechanism  is  the  dominant  factor  in 
Langmuir  circulations. 
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2.  Field  and  Laboratory  Observations  of  Langmuir  Circulations 

Basic  features  of  Langmuir  circulations  examined  in  either  field  or  laboratory 
studies  are  windrow  spacing,  vertical  velocities  below  the  convergence  and  divergence 
zones,  surface  velocities,  depth  of  penetration,  formation  times  after  wind  speed  reaches 
the  significant  velocity,  realignment  time  after  the  wind  speed  shift,  and  angle  between 
the  surface  streaks  and  the  wind  direction.  In  addition  to  these  observed  features  of  the 
circulations,  most  field  experiments  measure  wind  speed,  water  temperature,  and  water 
density  data  to  relate  the  observed  features  to  the  wind  forcing,  thermal  forcing  and 
density  structures  respectively.  The  threshold  wind  speed  to  initiate  the  circulations  is 
commonly  reported  as  3  ms'^  (Leibovich  1983).  Myer  (1971)  observed  that  circulations 
were  likely  to  occur  at  lower  wind  speeds  in  thermally  unstable  conditions  than  would  be 
required  in  stable  conditions. 

a.  Streak  Spacing 

Langmuir  (1938)  explained  the  difficulties  of  defining  the  distance 
between  surface  streaks  because  of  the  occurrence  of  numerous  smaller  streaks  between 
the  well-defined  larger  streaks.  The  larger  surface  vortices  contain  smaller  and  shallower 
vortices  in  the  same  way  large  waves  have  smaller  waves  upon  them  (Langmuir  1938). 
Langmuir  (1938)  reported  the  spacing  in  Lake  George  to  be  between  5  to  25  m  and 
changing  seasonally  depending  upon  the  thermocline  depth,  and  also  reported  the  spacing 
in  the  North  Atlantic  from  100-200  m.  Leibovich  (1983)  explained  that  the  row  spacing 
in  the  ocean  is  greater  than  in  lakes  and  ranges  from  two  to  hundreds  of  meters.  He  also 
cited  some  studies  such  as  Assaf  et  al  (1971),  Williams  (1965),  Ichiye  (1967),  Scott  et  al 
(1969),  Myer  (1969,  1971),  and  Harris  and  Lott  (1973),  which  showed  a  dynamical 
cascade  from  smaller  scales  towards  dominant  larger  scales. 

Langmuir  (1938)  and  the  most  of  the  following  studies  related  the  spacing 
(L)  to  the  penetration  depth  of  the  circulations  or  mixing  layer  depth  (D).  Langmuir 
(1938)  suggested  that  the  circulations  penetrate  down  to  the  bottom  of  the  thermocline, 
which  is  deeper  in  the  fall  than  in  the  spring  and  early  summer,  and  found  larger  spacing 
during  the  times  when  the  thermocline  was  deeper.  Assaf  et  al  (1971)  found  spacing  of 
the  largest  streak  scales  to  be  280  m,  with  a  mixed  layer  depth  of  200  m,  which  shows 
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that  LID  =  IA  and  Maratos  (1971)  reported  L !  D  =  1.66 .  Smith  (1992)  found  that  the 
spacing  was  2  to  3  times  the  mixed  layer  depth  while  it  varied  from  40-60  m  due  to  the 
baroclinic  tide  modulation  of  mixed  layer  depth.  Smith  (1992)  also  found  that  the 
circulations  deepened  the  mixed  layer. 

Faller  and  Caponi  (1978)  reported  that  L  /  H  =  4.8  [1-exp  (-0.5  Lw  /  H)] 
during  their  laboratory  experiments,  where  H  is  the  mean  depth  of  the  water  in  the  wind- 
wave  tank  or  estimated  mixed  layer  depth  for  oceanic  data,  L  is  the  spacing  between 
convergence  zones,  and  Lw  is  the  wavelength  of  the  surface  waves.  Their  laboratory 
results  and  several  sets  of  oceanic  data  fit  this  equation.  They  noted  that  the  circulations 
begin  within  seconds  after  arrival  of  the  waves  and  concluded  that  the  scale  of  the  surface 
waves  is  important  in  forcing  the  circulations. 

Uzaki  and  Matsunaga  (2000)  examined  the  Langmuir  circulations  in  a 
laboratory  wind- wave  tank.  The  depth  of  the  tank  was  varied  from  2  cm  to  16  cm.  They 
observed  that  the  spacing  depend  upon  the  water  depth  rather  than  the  wave  parameters 
and  wind  speed,  and  reported  the  spacing  2.5  times  the  water  depth  in  a  shallow  water 
region.  Faller  and  Woodcock  (1964)  related  the  mean  spacing  to  the  wind  speed  and 
found  that  L  =  (4.8  s)  x  Uw,  where  Uw  is  the  wind  speed  10  m  above  the  sea  surface. 
Maratos  (1971)  reported  a  similar  relationship  with  L  =  0.1m-l-(4.8s)x  Uw-  Leibovich 
(1983)  stated  that  the  number  of  samples  and  the  range  of  wind  speed  used  in  these 
experiments  were  small  and  these  relationships  may  not  be  accurate.  Scott  et  al  (1969) 
attempted  to  find  a  similar  relationship  but  could  not  find  a  statistically  significant 
correlation  between  spacing  and  wind  speed. 

b.  Orientation 

The  orientation  of  the  streaks  was  generally  observed  with  the  axis  nearly 
aligned  with  the  wind  direction.  Different  studies  reported  different  orientation  angles 
between  the  wind  and  the  streaks.  Leibovich  (1983)  noted  that  the  windrows  were  nearly 
aligned  with  wind  direction  with  a  maximum  deviation  of  20°.  Langmuir  (1938),  Walther 
(1967),  Kenny  (1977)  found  that  windrows  were  aligned  within  a  few  degrees  of  the 
wind  direction  in  their  experiments  in  the  lakes,  and  the  same  result  was  found  by 
Williams  (1965),  Welander  (1963),  and  Maratos  (1971)  in  the  ocean. 
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Faller  (1964)  looked  for  a  systematic  deviation  from  the  wind  direction 
and  analyzed  oceanic  data  from  the  Northern  Hemisphere,  finding  consistent  deflection 
of  13°  ±  2°  to  the  right  of  the  wind  direction. 


Leibovich  and  A.  Tandon  (1993)  found  that  the  small  windrow  angles 
were  due  to  weak  density  stratifications.  The  largest  angle  that  they  found  was  5.3 
degrees  with  an  equal  likelihood  of  being  on  either  side  of  the  wind  direction.  They 
reported  that  the  speed  of  the  lateral  drift  associated  with  the  angular  deviation  increases 
with  increasing  stratification  strength. 

c.  Formation  Times 

Pollard  (1977)  indicated  that  Langmuir  circulations  appear  within  a  few 
tens  of  minutes  after  the  onset  of  winds  having  speeds  larger  than  3  ms'\  The 
reorientation  time  of  the  circulations  after  the  wind  direction  shift  was  studied  by  several 
scientists.  Maratos  (1971)  observed  that  Langmuir  circulations  reoriented  with  the  new 
wind  direction  within  2  to  4  minutes  after  the  rapid  shift  of  wind  direction.  Langmuir 
(1938)  noted  that  the  reorientation  time  was  within  20  minutes  after  a  wind  direction  shift 
in  Atlantic  Ocean.  Welander  (1963)  observed  that  8-m  spaced  rows  reoriented  within  10 
minutes  after  the  wind  shift,  and  the  reoriented  rows  had  the  same  spacing  as  the  previous 
ones.  Stommel  (1951)  found  the  smallest  reorientation  time,  and  reported  that  small-scale 
windrows  on  ponds  reoriented  within  1  -2  minutes  after  the  wind  shift.  Faller  (1981) 
suggested  that  small  scale  Langmuir  circulations  have  more  rapid  responses  to  wind 
shifts  than  larger  scale  Langmuir  circulations. 

d.  Strength  and  Persistence 

Plueddemann  et  al  (1996)  related  Langmuir  circulations’  strength  to  wind 


speed  and  Stokes  drift  and  formulated  the  strength  as 


1/2 


where  u. 


is 


the  water  friction  velocity  and  is  the  Stokes  drift  velocity.  They  also  observed  that  the 
Langmuir  circulations  were  maintained  for  up  to  24  hours  after  the  wind  relaxation  and 
hypothesized  that  the  waves  were  responsible  for  this  persistence,  because  the  waves  had 
much  slower  relaxation  time  than  the  winds  during  their  observations. 
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e.  Surface  and  Subsurface  Current  Structure 

Surface  current  structures  associated  with  Langmuir  circulations  were 
determined  from  the  collecting  time  of  the  surface  tracers  into  rows,  and  the  vertical 
velocities  were  determined  from  the  sinking  motions  of  the  surface  particles  until  the 
1980s.  Smith  (2001)  noted  that  the  velocity  components  of  the  Langmuir  circulations  had 
started  to  be  estimated  using  the  Doppler  shift  of  the  backscattering  sound  signals  in  the 
late  1980s. 

Leibovich  (1983)  reported  that  the  along- wind  horizontal  velocity 
component  measured  at  streaks  was  larger  than  those  between  them.  Ichiye  (1967) 
reported  that  along-wind  velocities  at  streaks  were  two  times  larger  than  the  along-wind 
velocities  between  the  streaks.  10  cms'^  and  6  cms'^  differences  were  reported  by  Assaf  et 
al  (1971)  and  Harris  and  Lott  (1973)  respectively,  but  Leibovich  (1983)  stated  that  these 
results  were  reported  without  reference  to  other  components  of  the  current  system  or 
prevailing  wind  speed. 

Langmuir  (1938)  reported  average  cross  wind  horizontal  velocities  at  the 
surface  of  3  cms'^  while  Woodcock  (1944)  and  Ichiye  (1967)  reported  approximately 
6  cms'^  and  5  cms"'  respectively.  Plueddemann  et  al  (1996)  related  the  crosswind  surface 
velocities  to  the  wind  speed  and  waves,  and  parameterized  with  (WUs)  where  W  is 
wind  speed.  Smith  (1996)  parameterized  the  same  relation  with  (W  Us)  .  Both 
parameterizations  were  made  for  the  initial  growth  of  the  Langmuir  circulations. 
However  Smith  (2001)  examined  the  data  from  MBLEX-1  (marine  boundary  layer 
experiment)  and  noted  that  wind  speed  contribution  to  the  scaling  vanished  after 
Langmuir  circulations  had  developed,  and  described  the  scaling  as  v  ~  Us.  Smith  (2001) 
also  investigated  if  this  scaling  was  valid  for  other  data  sets  and  applied  this  scaling  to  the 
SWAPP  (Surface  Waves  Process  Program)  data.  He  found  that  the  relationship  between 
cross  wind  surface  velocities,  wind  speed,  and  Stokes  drift  was  still  statistically 
significant  but  correlations  were  not  as  strong  as  the  MBLEX-1  data. 

Eangmuir  (1938)  reported  downwelling  vertical  velocities  of  2  to  3  cms'^ 
below  the  surface  streaks,  and  upwelling  vertical  velocities  below  the  midway  adjacent 
streaks  were  1  to  1.5  cms'\  Myer  (1969,  1971)  found  the  downwelling  velocities  3  cms'^ 
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under  thermally  stable  conditions  and  5  cms"'  under  thermally  unstable  conditions. 
Leibovich  (1983)  described  maximum  downwelling  velocities  as  roughly  one  percent  of 
the  wind  speed  or  one  third  to  one  quarter  of  maximum  wind-induced  surface  currents. 

Sutcliffe  et  al  (1963)  correlated  the  downwelling  velocities  with  wind 

—2  1  1 

speed  and  found  that  w  =  0.85  x  10  W  or  1  cms'  per  1  ms'  wind,  where  W  is  wind 
speed  and  w  is  downwelling  velocity. 


Sutcliffe  et  al. 
Harris  and  Lott 

Fiiotov,  RjonzNn 
ond  Zoycev 


Figure  2.10  Measured  downwelling  speeds  below  streaks  as  a  function  of  wind  speed. 

The  open  squares  and  circles  correspond  to  surface  heating,  closed 
symbols  to  surface  cooling.  The  dashed  line  is  a  least  square  linear  fit  to 
all  of  the  Sutcliffe  et  al  (1963)  and  Harris  and  Lott  (1973)  data  taken 
together  while  the  solid  line  is  the  corresponding  fit  to  all  of  the  Filatov  et 
al  (1981)  data  (From  Leibovich  1983). 


Leibovich  (1983)  correlated  the  vertical  down  wind  velocities  with  the 
wind  speed  by  using  the  data  collected  by  Harris  and  Lott  (1973),  Filatov  et  al  (1981), 
and  Sutcliffe  (1963)  and  produced  the  correlation  plot  in  Figure  2.10  and  concluded  that 
downwelling  velocity  increases  linearly  with  increasing  wind  speed. 

Smith  et  al  (1987.)  found  that  the  vertical  velocities  range  from  15  to  30 
cms'^  in  15  ms''  winds  or  about  1  or  2  %  of  wind  speed. 
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III.  SITE  DESCRIPTION  AND  DATA  COLLECTION 


A  broadband  acoustic  Doppler  current  profiler  (BADCP)  was  used  to  measure 
velocity  profiles  across  the  water  column  during  the  Coupled  Boundary  Layers  Air  Sea 
Transfer  (CBLAST)  project.  Data  were  collected  in  2003  between  year  day  228  and  360. 
The  BADCP  was  cabled  to  the  Air  Sea  Interaction  Tower  (ASIT),  which  provided  power 
and  sent  a  continuous  data  stream  to  a  data  logging  computer  at  the  Martha’s  Vineyard 
Coastal  Observatory  (MVCO). 

A.  THE  COUPLED  BOUNDARY  LAYERS  AIR  SEA  TRANSFER  (CBLAST) 

The  coupled  boundary  layers  air-Sea  transfer  (CBLAST)  project  focuses  on 
processes  that  occur  in  the  oceanic  and  atmospheric  surface  boundary  layers,  which  are 
regions  influenced  by  ocean  surface  waves.  This  research  was  sponsored  by  the  Office  of 
Naval  Research  (ONR)  to  develop  improved  parameterizations  to  be  used  in  high- 
resolution  coupled  air-ocean  regional  numerical  models. 
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Figure  3.1  The  Physical  Processes  in  the  atmospheric  and  oceanic  boundary  layers 
examined  in  CBLAST  (From  www.whoi.edu) 


35 


CBLAST  investigated  the  processes  that  couple  momentum,  mass,  and  heat  across  the 
turbulent  atmospheric  and  oceanic  boundary  layers  in  the  presence  of  surface  gravity 
waves.  The  main  objectives  of  CBLAST  were  to  understand  the  influence  of  surface 
waves  in  air-ocean  coupling  and  determining  the  3-D  structures  of  these  boundary  layers, 
and  to  improve  the  understanding  of  air-sea  interactions  in  low  to  moderate  winds. 

Figure  3.1  summarizes  processes  studied  in  the  CBLAST  project.  CBLAST  had 
both  observational  and  modeling  components.  The  observational  measurement  tools  are 
fixed  towers  and  moorings,  satellite  and  airborne  remote  sensing,  autonomous  under 
water  vehicle  (AUV),  research  vessels  and  aircraft.  The  in  situ  observations  included 
oceanic  and  atmospheric  turbulence  and  mean  flow  from  the  fixed  towers  and  moorings, 
2-D  and  3-D  structure  of  the  boundary  layers,  2D  ocean  surface  structure  from  the  remote 
sensing  measurements,  and  surface  wave  measurements  with  a  special  emphasis  on  small 
scale  and  wave  breaking.  All  the  observational  investigations  will  help  quantify  the 
turbulent  kinetic  energy  (TKE)  momentum,  mass  and  heat  budgets  in  the  oceanic  mixed 
layer  and  atmospheric  boundary  layer  in  a  hierarchy  of  large  eddy  simulation  (LES)  and 
regional  coupled  atmosphere  /  ocean  models. 

1.  Air  Sea  Interaction  Tower  (ASIT) 

The  ASIT  (Figure  3.2)  was  built  to  meet  some  of  the  observational  requirements 
for  CBLAST.  ASIT  was  instrumented  from  top  to  bottom  with  sensors  capable  of 
measuring  momentum,  heat  and  mass  exchange  between  the  atmosphere  and  the  ocean. 
The  tower  structure  was  designed  by  Appledore  Engineering  Inc.  and  constructed  by 
RDA  Construction  Corporation  under  the  direction  of  Dr.  Tim  Edson  at  WHOI,  and  was 
funded  by  ONR.  It  was  specifically  designed  to  make  the  measurements  in  the  ocean  and 
the  atmosphere  to  investigate  the  ocean  processes,  which  include  air-sea  interactions,  gas 
exchange,  ocean  mixing,  bio-optics  and  sediment  transport. 
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Figure  3.2  Air  Sea  Interaction  Tower  (Fromwww.whoi.edu) 

The  ASIT  is  located  3.5  km.  due  south  of  MVCO  (Figure  3.3)  in  16.5  m.  water 
and  extends  25  m.  above  the  water  line  into  the  marine  boundary  layer.  It  is  connected  to 
Martha’s  Vineyard  coastal  observatory  (MVCO)  by  a  hybrid  fiber  optic  data  and  power 
cable.  Routine  data  measured  by  ASIT  are  integrated  with  MVCO  data,  and  the 
integrated  data  time  series  includes  wind  speed  and  wind  direction,  tidal  heights,  air  and 
sea  temperature,  off  shore  bottom  and  surface  currents,  and  wave  directional  spectra.  The 
ASIT  location  is  41°  19.5'  N  70°  34.0'  W  and  marked  with  the  letter  C  in  Figure  3.3. 


A:  MVCO  Meteorological  Mast 
B:  MVCO  Offshore  Node 
C:  Air-Sea  Interaction  Tower 


Figure  3.3  Martha’s  Vineyard  with  gray  color  and  ASIT  Location  marked  with  ‘C’ 
symbol  3.5  km  south  of  Martha’s  Vineyard  (From  www.whoi.edu) 
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The  BADCP  used  in  this  study  was  deployed  by  Res.  Prof.  T.  P.  Stanton  as  a 
component  of  his  ocean  measurements  for  CBLAST. 


Figure  3.4  BADCP  and  ASIT  orientation  at  the  ocean  bottom.  The  distance  between 
ASIT  and  the  BADCP  is  about  20m  and  there  is  data  node  between  the 
ASIT  and  BADCP 

The  NPS  5-beam  BADCP  was  deployed  offshore  20  meters  south  of  ASIT  from  a  pipe 
jetted  into  the  seabed,  looking  upward  to  the  water  surface  (Figure  3.4).  The  distance 
between  the  BADCP  and  the  seabed  was  about  1.5  m  and  the  first  bin  was  1  m  above  the 
transducer.  The  BADCP  was  connected  to  an  underwater  data  node,  which  provided 
communication  to  the  shore-based  Linux  data-logging  computer  through  the  ASIT 
ethemet  system.  The  -0.5  Hz  sampled  profile  data  was  acquired  by  a  C-program  running 
on  the  shore  computer  at  the  MVCO  shore  lab,  and  the  data  files  were  routinely 
transferred  back  to  NPS  for  processing  and  analysis.  The  local  data  logging  removed 
dependence  on  the  T1  data  link  from  MVCO  to  the  mainland. 

B.  BROADBAND  ACOUSTIC  DOPPLER  CURRENT  PROFILER  (BADCP) 

The  custom-built  5-beam  BADCP  manufactured  by  RD  Instruments  transmits 
acoustic  pulses  from  five  individual  transducers  along  five  individual  beams,  and  the 
transducers  receive  the  sound  backscattered  from  the  zooplankton  and  small  particles 
such  as  bubbles  floating  with  the  water  currents.  By  comparing  the  frequency  of 
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transmitted  and  received  sound  signals,  the  BADCP  measures  the  velocity  profiles  along 
the  beams  at  each  depth  cell  by  using  the  Doppler  shift  principle. 


Figure  3.5  Typical  BADCP  profiling  technique  that  allows  collecting  the  data  at 
different  depth  cells  (bins)  along  five  beams.  Depth  cells  represent 
conventional  single-point  current  meters  (From  RDI  BADCP  technical 
manual) 

The  depth  cell  shown  in  Figure  3.5  is  the  vertical  length  of  each  cell  or  bin  where  the 
velocities  are  measured.  Depth  cells  or  bins  represent  the  conventional  single  point 
current  meters  that  are  deployed  in  different  depths  throughout  the  water  column.  The 
BADCP  not  only  collects  the  velocity  data,  but  also  calculates  the  echo  intensity  of  the 
backscattered  sound  signals  and  the  correlation  magnitudes  at  each  bin.  Echo  intensity  is 
useful  in  determining  signal  strength  and  the  relative  volume  of  the  particles 
backsc altering  the  sound  signal  to  the  transducers,  and  can  be  used  to  check  data  quality 
too.  Correlation  magnitudes  are  the  measure  of  the  pulse-to-pulse  correlation  in  a  ping  for 
each  bin,  and  represent  the  validity  or  confidence  of  the  data. 

Broadband  ADCPs  transmit  a  series  of  short  pulses  instead  of  one  long  pulse. 
Since  this  feature  allows  the  separation  of  pulse  width  and  pulse  lag,  longer  lags  can  be 
chosen  for  lower  velocity  noise  measurements  independent  of  the  size  of  the  pulse. 
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Figure  3.6  Phase  change  process  for  the  Broadband  System  (From  RDI  BADCP 
technical  manual) 

In  Figure  3.6,  a  broadband  ADCP  transmits  two  pulses  of  duration  Tp,  separated  by  a  lag 
Tp.  Figure  3.6-A  shows  no  phase  change,  Figure  3.6-B  shows  a  small  phase  change  and 
Figure  3.6-C  shows  a  bigger  phase  change  due  to  the  relative  velocity  of  the  particles, 
which  backscatter  the  transmitted  sound  signal  to  the  transducers. 


Figure  3.7  NPS  five-beam  BADCP  Beam  geometry  (From  RDI  BADCP  technical 
manual) 

The  BADCP  coordinate  system  and  the  transducer  orientation  is  shown  in  Figure 
3.7-A.  Beam  1  and  beam  2  measure  the  combination  of  the  velocity  components  in  the  x 
and  z  directions,  while  beam  3  and  beam  4  measure  the  combination  of  the  velocity 
components  in  the  y  and  z  directions.  Beam  5  only  responds  to  vertical  velocity 
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components.  Figure  3.7-B  is  the  side  view  of  any  two  slanted  beams  and  shows  measured 
components  of  the  aetual  eurrent. 


The  BADCP  has  four  slanted  beams  and  one  vertical  beam  and  the  angle  between 
the  slanted  beams  and  the  vertical  is  30  degrees. 


BADCP 

Depth  Cell  or 

Number  of 

Pings  Per 

Ensemble 

Frequeney 

Bin  Size 

Depth  Cells 

Ensemble 

Interval 

600  kHz. 

0.5  m. 

40 

10 

-  2  sec. 

Table  3.  NFS  five  BADCP  parameters  used  during  CBLAST  experiment. 


The  instrument  parameters  used  for  this  experiment  are  shown  in  Table  3.  The  ensemble 
interval  depends  on  the  pings  per  ensemble.  When  pings  per  ensemble  are  adjusted  to  10, 
the  ensemble  interval  is  automatically  assigned  by  the  BADCP  to  be  about  two  seconds 
providing  a  -'0.5  Hz  data  rate. 
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IV.  DATA  PROCESSING 


A  series  of  matlab  program  modules  were  written  to  load  the  0.5  Hz  (two-second) 
beam  profile  data,  find  the  surface  return  bin,  and  determine  thirty-second  low-pass 
filtered  profile  time  series  of  velocity  and  backscatter  data  beneath  the  surface.  An 
algorithm  to  detect  the  upwelling  and  downwelling  regions  of  Langmuir  circulations  was 
created  based  on  the  coherency  features  of  these  circulations  using  the  thirty-second 
beam  five  vertical  velocities  in  two-hour  intervals.  Following  cell  detection,  maximum 
upwelling  and  downwelling  velocities,  penetration  depth  and  spacing  were  defined. 

A.  LOADING  AND  BLANKING  THE  DATA 

1.  Loading  the  0.5  Hz  Data 

The  Loadbadcp  matlab  function  was  created  to  load  the  two-second  data. 
Loadbadcp  uses  the  mex  functions  findfile  and  arin,  which  link  into  the  raw  database. 
Findfile  identifies  files  for  each  desired  time  period,  and  arin  loads  the  velocity,  intensity, 
and  correlation  data  from  each  of  these  files.  Loadbadcp  concatenates  the  loaded  data  into 
the  variables  assigned  for  each  data  type.  For  example,  vl,  v2,  v3,  v4  and  v5  are  the 
variables  for  the  velocity  data  of  beam  1,  beam  2,  beam  3,  beam  4  and  beam  5 
respectively.  Loadbadcp  also  used  the  pdrin  matlab  function  to  get  the  header  information 
of  the  BADCP  data,  including  tilt  and  heading  information  for  the  deployed  instrument. 
In  the  following  analysis  thirty-second  segments  of  0.5  Hz  data  were  loaded. 

2.  Blanking  the  Data 

The  BADCP  used  in  this  experiment  has  no  pressure  sensor  in  it,  making  it 
difficult  to  define  the  sea  surface  accurately,  since  the  sea  surface  is  changing 
continuously  in  response  to  both  waves  and  tides.  However  since  the  BADCP  was 
programmed  to  collect  the  data  for  forty  bins  (20  m),  well  beyond  the  16  m  instrument 
depth,  the  sea  surface  was  determined  from  the  strong  surface  returns  seen  in  the  intensity 
profile  time  series  of  beam  5.  This  allowed  the  reference  bin  number,  or  top  safe  bin 
closest  to  the  sea  surface,  to  be  defined  for  each  0.5  Hz  raw  profile  sample. 

After  loading  the  two-second  data,  the  sea  surface,  corresponding  to  the  bin 
number  of  maximum  intensity,  was  found  between  bins  20  and  40  (Figure  4.1).  The 
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vertical  beam  5  was  chosen  for  this  process  to  avoid  the  side  lobe  and  penetrating  bubble 
effects  seen  in  the  sloping  beams. 


Figure  4.1 


The  vertical  beam  5  intensity  time  series  color  plot  for  40  seconds  during 
year  day  318 


Once  the  surface  bin  was  selected,  the  bin  number  three  below  it  was  chosen  as  the 
blanking  boundary,  and  all  the  data  above  this  blanking  boundary  was  replaced  by  NaNs 
(NaN  means  not  a  number  in  the  Matlab  programming  language).  This  process  was  done 
for  each  ensemble  of  two-second  data  for  thirty-second  time  periods  throughout  the  data 
set.  Vertical  beam  5  data  were  directly  blanked  while  slanted  beams  were  not. 


B.  HORIZONTAL  AND  VERTICAL  VELOCITY  COMPONENTS 

Horizontal  velocities  were  calculated  with  two  different  methods,  a  four-beam 
solution  of  slanted  beams  and  a  two-beam  solution  exploiting  the  vertical  beam  5. 

1,  Four-Beam  Solution  of  Horizontal  Velocity  Components 

The  four-beam  solution  of  horizontal  velocity  components  is  typically  used  in 
four-beam  ADCP’s.  The  horizontal  velocity  component  in  the  x-direction,  u,  is  calculated 
by  differencing  beam  1  and  beam  2  (Figure  4.2-A),  while  the  horizontal  velocity 
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component  in  the  y-direction,  v,  is  ealeulated  by  differeneing  beam  3  and  beam  4  (Figure 
4.2-B).  All  these  solutions  presume  that  the  horizontal  flow  is  uniform  over  the  spatial 
seale  between  the  seleeted  beams  at  each  bin  depth. 


Figure  4.2-A  Figure  4.2-B 

Figure  4.2  Horizontal  veloeity  eomponents  four-beam  solution  geometry.  The  u- 
eomponent  is  calculated  between  beam  1  and  beam  2  (Panel  A),  the  v- 
component  is  calculated  between  beam  3  and  beam  4  (Panel  B)  (From 
RDI  BADCP  technical  manual) 

The  formulations  for  the  four-beam  solution  are  based  on  those  found  in  RDFs 
publication,  “  Direct  Reading  and  Self-Contained  Broadband  Acoustic  Doppler  Current 
Profiler  Technical  Manual”. 


u  =  (-Fl  +  F2)/sin6’ 

(4.1) 

v  =  (-F3  +  F4)/sin6' 

(4.2) 

where,  VI,  V2,  V3,  V4  are  the  measured  beam  1,  beam  2,  beam  3,  beam  4  velocities 
respectively  and  0  is  the  thirty  degree  angle  between  the  vertical  axis  and  the  slanted 
beams.  These  formulas  assume  a  uniform  flow  across  the  aperture  formed  by  the 
diverging  beams  at  each  bin. 
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2,  Two-Beam  Solution  of  Horizontal  Velocity  Components 

The  two-beam  solution  method  calculates  the  horizontal  velocity  components 
between  the  slanted  beams  and  vertical  beam  by  using  the  trigonometry  between  beams 
as  shown  earlier  in  Figure  3.7-B. 


Figure  4.3-A  Figure  4.3-B 

Figure  4.3  Panel  A;  Two  different  calculation  of  the  u-component  with  the  two-beam 
solution  between  beams  1  and  5,  and  beams  2  and  5.  Panel  B;  Two 
different  calculations  of  the  v-component  with  the  two-beam  solution 
between  beams  3  and  5,  and  beams  4  and  5.  The  lateral  scale  over  which 
the  horizontal  velocity  estimates  are  made  is  1/2  that  of  the  four  beam 
solution. 

The  two-beam  solution  (Figure  4.3)  allows  higher  resolution  estimates  of  u  and  v 
than  the  four-beam  solution  as  it  uses  half  the  scale  of  the  four-beam  solution.  The  u- 
component  between  beam  1  and  beam  5  is  given  by 

u  =  -[[V,-{V,x  Cose)) /Sine)  (4.3) 

while  the  u-component  between  beam  2  and  beam  5  is  given  by 

M  =  (jz  - (F5  X  Cose)) / Sine  (4.4) 

and  the  v-component  between  beam  3  and  beam  5  is  given  by 

v  =  -[[V,-{V,x  Cose))/ Sine)  (4.5) 
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while  the  v-eomponent  between  beam  4  and  beam  5  is  given  by 

v  =  (V,-(V,xCos0))/Sm0  (4.6) 

where  Vi,  V2,  V3,  V4  and  V5  represent  the  measured  beam  1,  beam  2,  beam  3,  beam  4  and 
beam  5  velocities  respectively,  and  ^  is  the  thirty  degree  angle  between  the  slanted 
beams  and  vertical  beam. 

3.  Vertical  Velocity  Solutions 

The  fifth  vertical  beam  produces  vertical  velocity  profiles  from  -0.2  m  diameter 
by  0.5  m  vertical  sample  volumes,  and  consequently  does  not  require  the  horizontal 
homogeneity  presumption  seen  in  the  u  /v  velocity  estimates  above.  The  vertical  fifth 
beam  also  provides  two  u,  and  two  v  estimates  as  discussed  above,  allowing  crude 
horizontal  strain  estimates  to  be  made  at  bin-lag  scales. 

C.  TIME  AVERAGING  THE  DATA 

After  blanking  the  two-second  beam  5  vertical  velocities  and  calculating 
horizontal  velocity  components,  the  data  were  averaged  into  thirty-second  blocks  to  reject 
orbital  motion  from  surface  gravity  waves  using  the  matlab  function  zapmeanmat. 
Zapmeanmat  is  an  NFS  turbulence  group  function  that  calculates  the  mean  of  every  row 
in  an  M  by  N  matrix  ignoring  NaNs,  which  were  put  in  the  velocity  matrices  by  the 
blanking  process.  (Figure  4.4) 


Figure  4.4 


NaN 

NaN 

NaN 

NaN 

NaN 

i4 

NaN 

NaN 

6 

3 

NaN 

6 

4 

2 

1 

3 

3 

3 

1 

2 

2 

/..\i>Mi;.\NM.vr 


NaN 

4 

5 
3 
2 


Thirty-second  temporal  averaging  of  the  BADCP  data.  The  5  by  4  matrix 
represents  the  40  by  15  raw  data  matrix.  The  rows  are  numbered  from 
bottom  to  top  and  represent  the  bin  numbers  while  the  columns  represent 
the  two-second  ensembles,  with  15  columns  of  data  loaded  for  every 
thirty-second  time  period. 


Loading,  blanking  and  averaging  processes  were  done  in  the  “mnvelcom”  matlab 


routine.  Mnvelcom  loops  through  two-hour  intervals  between  selected  start  and  end  days. 
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It  loads  the  raw  data  in  thirty-second  increments  at  each  step  of  the  loop,  and  applies  the 
blanking  and  averaging  processes.  Finally,  mnvelcom  creates  thirty-second  beam  5 
vertical  velocities,  horizontal  velocity  components  and  intensity  profile  time  series  for 
each  two-hour  period.  Mnvelcom  also  calculates  the  mean  and  the  standard  deviation  of 
the  boundary  bin  numbers  of  the  blanking  process  for  every  thirty-second  period  in  the 
loop,  as  shown  in  Figure  4.5 


Figure  4.5  Mnvelcom  matlab  routine  algorithm  for  loading  the  0.5  Hz  data,  blanking 
and  thirty-second  temporal  averaging  processes  and  creating  thirty-second 
data  for  two-hour  periods. 

D.  DEFINING  THE  NEAR-SURFACE  REFERENCE  BIN  NUMBER 

A  near-surface  reference  bin  number  was  found  for  the  bin  closest  to  the  sea 
surface  that  has  valid  data  for  each  thirty-second  data  set  in  every  two-hour  period.  The 
Refbindefin.m  matlab  routine  was  created  to  define  the  reference  bin  number,  which 
allows  near-surface  vertical  velocity  to  be  calculated  for  detecting  Langmuir  circulations 
in  subsequent  analyses. 

The  near-surface  reference  bin  number  was  found  from  profiles  of  standard 
deviations  of  the  thirty-second  averaged  data  for  each  two  hour  period.  Under  energetic 
waves  with  large  amplitudes  (Figure  4.6),  there  is  a  distinct  break  in  vertical  velocity 
standard  deviations  increasing  toward  the  surface  due  to  surface  beat,  turbulence  and 
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Langmuir  circulations  and  very  high  values  close  to  the  surface  (Figure  4.7).  These  large 
standard  deviation  values  are  attributed  to  surface-hit-contaminated  vertical  velocity 
measurements  not  screened  out  from  the  relatively  slow  two-second  sampling  of  intensity 
profiles.  This  distinct  breakpoint  allowed  these  velocity  bins  to  be  excluded  from 
subsequent  analyses. 


Figure  4.6  Wind  speed  and  significant  wave  heights  (SWFI)  for  the  days  between  350 
and  354.  The  energetic  high  amplitude  waves  are  in  between  days  351.5 
and  353.5 
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Figure  4.7  Temporal  Standard  Deviations  (dots)  of  beam  5  vertical  velocities  for  the 
high  wind  /  wave  time  period  between  the  days  352  and  352.0833.  Solid 
line  is  the  mean  bin  number  (or  mean  sea  level)  defined  during  the 
blanking  process  while  red  and  green  lines  are  standard  deviations  from 
mean  sea  level. 

To  illustrate  this  screening  process,  Figure  4.6  shows  the  significant  wave  heights  and 
wind  speed  between  the  days  350  and  354  with  wave  height  increasing  from  0.8m  to 
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3.5m  during  twelve-hour  period.  Figure  4.7  shows  the  standard  deviations  (dots)  of  the 
thirty-second  beam  5  vertical  velocities  for  a  two-hour  period  between  year  day  352  and 
352.0833. 

The  slope  between  each  dot  was  calculated  as 

,  A{std(V5)) 

slope  =  — ^ i — —  (4.7) 

A(Bin  No) 

where  std  is  the  standard  deviation  in  time,  V5  is  beam  5  vertical  velocities  in  ms'^  and  A 
represents  the  difference  between  each  point  and  each  bin.  For  the  slopes  greater  than 
0.01,  the  corresponding  velocity  bins  were  considered  as  bad.  The  bin  number  one  below 
the  highest  bin  whose  slope  was  less  than  0.01  was  chosen  as  reference  bin  number  for 
high  wave  cases.  For  the  case  shown  in  Figure  4.7,  the  reference  bin  number  was  defined 
as  19. 

For  low  wave  cases  (Figure  4.8),  the  same  method  was  used  as  the  high  wave 
case  but  in  this  case  the  bin  number  two  below  the  slope  break-point  was  chosen  as 
reference  bin  number. 


Time  (year  day) 


Figure  4.8  Wind  speed  and  significant  wave  heights  between  (SWH)  the  days  328 
and  332  representing  the  low  wave  case. 
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Figure  4.9  Temporal  Standard  Deviations  (dots)  of  beam  5  vertical  velocities  for  low 
wind  /  wave  time  period  between  the  days  330  and  330.0833.  Solid  line  is 
the  mean  bin  number  (or  mean  sea  level)  defined  during  the  blanking 
process  while  red  and  green  lines  are  standard  deviations  from  mean  sea 
level. 

Figure  4.9  represents  the  standard  deviations  of  the  beam  5  velocities  for  the  time  period 
between  the  year  day  330  and  330.0833  to  show  a  typical  low  wind  and  low  wave  case. 
Low  winds  and  low  waves  are  dominant  in  day  332  (Figure  4.8).  If  the  standard  deviation 
of  the  blanking  boundary  bin  numbers  shown  with  green  or  red  dashed  lines  from  mean 
sea  level  (solid  line)  in  Figure  4.9  is  less  than  0.47  (about  0.2  m),  the  program  considers 
this  case  as  low  wave  case  and  defines  a  reference  bin  number  two  bins  below  the  break¬ 
point  bin  number  estimated  by  the  algorithm  explained  above. 

Having  defined  the  near-surface  top  safe  velocity  bin,  all  the  data  matrices’  sizes 
were  adjusted  to  their  reference  bin  number.  For  example,  if  the  reference  bin  number 
was  20,  bin  numbers  above  the  twentieth  bin  were  neglected  in  a  40  by  240  matrix. 
Therefore  adjusted  the  matrix  became  20  by  240. 

E.  FILTERING  THE  DATA 

The  beam  5  vertical  velocities  and  the  horizontal  velocity  components  were  two 
minute  low-pass  filtered  in  time  to  minimize  the  residual  effects  of  orbital  swell  motion. 
The  vertical  velocity  data  were  also  20-minute  high-pass  filtered  in  time  to  reject  the 
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horizontal  velocities  entering  from  the  slightly  tilted  vertical  beam.  The  NFS  Lpfilt 
matlab  function  was  used  to  apply  low-pass  filtering  while  hpfilt  matlab  function  was 
used  to  apply  high-pass  filtering,  with  both  filters  using  a  four-pole  Butterworth  response. 
Beam  5  vertical  velocity  data  also  was  spatially  low-pass  filtered  vertically  through  the 
water  column  with  a  4  m  cutoff  and  second  order  filter,  before  calculating  the  maximum 
upwelling  /  downwelling  velocities  and  the  penetration  depths  of  the  Langmuir 
circulations. 

F.  LANGMUIR  CIRCULATIONS  DETECTION  PROCESS 

After  completing  the  processes  described  above,  beam  5  vertical  velocities  were 
used  to  detect  the  Langmuir  circulations.  The  detection  process  was  based  on  the 
coherency  features  of  the  Langmuir  circulations  through  the  water  column. 

The  mean  of  the  top  four  clean  velocity  bins  was  calculated,  where  the  top  of 
these  four  bins  was  the  reference  bin  number.  The  near-surface  mean  velocities  were  used 
to  define  upwelling  and  downwelling  regions  of  Langmuir  circulations. 
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Figure4.10  Langmuir  circulations  Detection  Process:  The  color  contour  plot  of 
vertical  velocities  (Panel  A);  Time  series  plot  of  upper  four  bin  (2  m) 
vertical  velocities  where  the  reference  bin  number  is  19  (Panel  B);  Mean 
of  the  upper  four  bin  (2m)  vertical  velocities  and  finding  all  the  turning 
points  where  the  circles  are  downwelling  velocities,  crosses  are  upwelling 
velocities  (Panel  C);  Applying  thresholds  +  0.005  ms’'  for  upwelling 
velocities,  and  -0.01  ms’'  for  downwelling  velocities  in  panel  C  and  yields 
panel  D 
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Figure  4.10  represents  the  beam  5  vertieal  veloeity  profile  ineluding  the  upwelling  and 
downwelling  veloeities  between  the  days  352  and  352.0833  (two  hour  period)  when  the 
high  waves  and  high  winds  were  dominant.  In  this  ease,  referenee  bin  number  is  19. 
Figure  4.10-A  is  a  eolor  eontour  plot  of  beam  5  velocities  after  the  band  pass  filtering 
process,  with  the  color  scale  on  the  right  panel.  Downwelling  regions  are  presented  by  the 
bluish  colors  while  the  upwelling  regions  presented  by  reddish  colors.  Figure  4.10-B  is 
the  top  four  bins’  velocities  of  beam  5  versus  time.  Each  bin  is  represented  by  a  different 
color.  In  Figure  4.10-C,  the  solid  blue  line  is  the  mean  of  the  upper  two-meter  clean 
vertical  velocities  shown  in  Figure  4.10-B.  The  plus  signs  in  Figure  4.10-C  show  the 
upward  turning  points  of  the  solid  blue  line  and  the  circles  show  the  downward  turning 
points  of  the  solid  blue  line,  and  thus  show  the  upwelling  and  downwelling  velocities 
respectively.  In  order  to  detect  the  Langmuir  circulations  velocity  thresholds  were  applied 
to  these  upwelling  and  downwelling  regions.  Velocities  less  than  +  0.005  ms'^  were 
neglected  for  upwelling  regions  while  the  velocities  less  than  -  0.01  ms'^  were  neglected 
for  downwelling  regions.  The  plus  signs  and  the  circles  in  Figure  4.10-D  show  the 
remaining  upwelling  and  downwelling  regions  after  velocity  thresholds  respectively. 

A  vertical  coherency  test  was  applied  to  the  remaining  points  as  the  last  step  of 
the  detection  process.  The  vertical  velocities  of  the  top  four  bins  were  tested  for 
coherency  by  measuring  the  normalized  standard  deviation  (NSD)  of  the  top  four  clean 
bins  at  each  selected  upwelling  /  downwelling  event.  For  NSD  >  0.35  across  these  four 
bins,  the  detected  event  was  deemed  vertically  incoherent,  and  was  rejected.  The  value  of 
0.35  was  chosen  after  many  trials  of  data  sets  for  different  cases,  including  low  winds- 
low  waves,  high  winds-high  waves,  increasing  winds-increasing  waves,  and  decreasing 
winds-  decreasing  waves.  For  NSD  less  than  0.35  corresponding  downwelling  and 
upwelling  points  were  assigned  as  Langmuir  circulation  events. 
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Normalized  Standard  Deviation  Distance  Below  Refbin  (m) 


Normalized  Beam  5  velocities  for  Downwelling  Regions 


-A- 


Normalized  Standard  Deviations  of  Top  Four  Bin  Velocities  for  Downwelling  Regions 


-B- 

Figure  4.1 1  Normalized  Beam  5  vertical  velocities  (Panel  A)  and  Normalized  standard 
deviations  (NSD)  of  top  four  bin  vertical  velocities  (Panel  B)  for 
downwelling  events  shown  in  Figure  4.10-D.  The  horizontal  solid  line  in 
Panel  B  represents  the  0.35  NSD  thresholds. 
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Normalized  Standard  Deviation  Distance  Below  Refbin  (m) 


Normalized  Standard  Deviations  of  Top  Four  Bin  Velocities  for  Upwelling  Regions 
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Figure  4.12  Normalized  Beam  5  vertical  velocities  (Panel  A)  and  Normalized  standard 
deviations  (NSD)  of  top  four  bin  vertical  velocities  (Panel  B)  for 
upwelling  events  shown  in  Figure  4.10-D.  The  horizontal  solid  line  in 
Panel  B  represents  the  0.35  NSD  thresholds. 
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In  Figures  4.1 1  and  4.12  the  top  four  bin  normalized  veloeities  for  eaeh  downwelling  and 
upwelling  points  have  the  similar  eolor  shading  if  their  NSD  values  are  less  than  0.35 
while  they  have  distinguishable  eontrast  eolors  (i.e.  red,  yellow,  blue  etc),  if  their 
correspondent  NSD  values  are  greater  than  0.35. 
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Filtered  Thirty-second  Beam  5  velocity  (m  s  ) 


Top  Four  bins  velocities  of  Beam  5 


Top  Four  Bins  Mean  Velocities  of  Beam  5  With  Turning  Points 
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The  Upwelling  and  Downwelling  Points  after  Velocity  Threshold 


The  Upwelling  and  Downwelling  Points  After  Coherency  Test 


Figure  4.13  Langmuir  circulations  detection  process:  Detected  downwelling  and 
upwelling  velocities  of  Langmuir  circulations  after  the  coherency  test  of 
the  downwelling  and  upwelling  regions  shown  in  Figure  4.10  -E.  The  first 
four  panels  are  the  same  as  in  Figure  4.10. 
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The  detection  algorithm  was  applied  to  other  cases  with  different  forcing  conditions  as 
follows. 


Forcing  Mechanisms  and  Ocean  Response  with  BADCP  Data  Time  Lines 
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Figure  4.14  Atmospheric  and  oceanic  forcing  mechanisms  with  the  available  BADCP 
data  in  winter  time  between  the  days  326  and  334.  The  wind  speed  in 
Panel  A  is  low  to  moderate  till  the  day  332.3  and  after  this  day  storm 
system  is  effective.  The  increasing  wind  speed  from  0  to  14  ms'^  between 
the  days  332.3  and  333.5,  trigger  the  ocean  and  wave  heights  in  Panel  B 
increases  from  0  to  4  m  during  the  same  time  period  which  is  a  good  case 
for  Langmuir  circulation  development.  The  wave  period  is  in  Panel  C  and 
the  tidal  currents  at  the  surface  in  Panel  D,  which  are  strong  enough  to 
compete  with  the  Langmuir  circulations  horizontal  velocities.  Times  of 
available  BADCP  data  are  in  Panel  E  for  this  time  period  of  year. 
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Two  hour  time  periods  between  the  days  326  and  334  (Figure  4.14)  were  chosen  to  apply 
the  detection  algorithm  for  different  forcing  mechanisms.  The  following  5  figures 
represent  the  detected  downwelling  /  upwelling  regions  of  Langmuir  circulations  and 
have  the  same  format  of  Figure  4.13. 
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The  Upwelling  and  Downwelling  Points  after  Velocity  Threshold 
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The  Upwelling  and  Downwelling  Points  After  Coherency  Test 


At  A: 


A 


328.01  328.02  328.03  328.04  328.05  328.06  328.07  328.08 

Time  (year  day) 

-E- 

Figure  4.15  Detection  of  upwelling  /  downwelling  regions  of  Langmuir  circulations 
between  the  days  328  and  328.0833.  The  wind  speed  is  ~2  ms'^  and  the 
wave  height  is  ~  0.3  m.  6  upwelling  regions  (Panel  E)  were  detected. 
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The  vertical  velocities  (Figure  4.15-A)  are  low  as  well  as  the  wind  speed  and  wave 
heights.  So  most  of  the  upwelling  /  downwelling  velocities  detected  in  Figure  4.15-C  is 
eliminated  after  the  velocity  thresholds  and  remaining  ones  are  in  Figure  4.15-D.  Only 
six  upwelling  regions  pass  the  coherency  test  (Figure  4.15-E)  for  this  low  wind  speed  and 
wave  height  case. 
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Figure  4.16  Detection  of  upwelling  /  downwelling  regions  of  Langmuir  circulations 
between  the  days  331  and  331.0833.  The  wind  speed  is  ~5  ms’'  and  the 
wave  height  is  ~  1  m.  12  upwelling  regions  and  11  downwelling  regions 
(Panel  E)  were  detected. 
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The  increasing  wind  speed  and  wave  height  case  (locally  developed  seas)  is  between  the 
days  331  and  331.0833  (Figure  4.16).  The  vertical  velocities  (Figure  4.16-A)  are  stronger 
and  there  are  more  coherent  upwelling  /  downwelling  velocities  through  the  water 
column  than  the  previous  case.  This  is  a  typical  case  for  Langmuir  circulations  to  occur  in 
the  ocean,  which  has  increasing  wave  height  due  to  increasing  wind  speed. 


Top  Four  bins  velocities  of  Beam  5 


The  Upwelling  and  Downwelling  Points  after  Velocity  Threshold 
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Figure  4.17  Detection  of  upwelling  /  downwelling  regions  of  Langmuir  circulations 
between  the  days  332  and  332.0833.  The  wind  speed  is  ~1  ms'^  and  the 
wave  height  is  ~  0  m.  5  upwelling  regions  (Panel  E)  were  detected.  This  is 
the  time  period  right  before  a  storm  case  between  the  days  332.6  and  334. 
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The  case  shown  in  Figure  4.17  has  low  wind  speed  /  wave  height  and  few  detected 
upwelling  velocities  (Figure  4.17  -E).  This  is  a  time  period  right  before  a  storm  case.  The 
following  two  cases  are  for  developing  and  well-developed  storm  cases  to  show  the 
increasing  amount  of  detected  Langmuir  circulations. 
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Figure  4.18  Detection  of  upwelling  /  downwelling  regions  of  Langmuir  circulations 
between  the  days  333  and  333.0833.  The  wind  speed  is  ~7  ms’'  and  the 
wave  height  is  ~  1.2  m.  20  upwelling  regions  and  9  downwelling  regions 
(Panel  E)  were  detected.  The  beginning  stage  of  the  storm  event  where  the 
wind  speed  and  wave  height  are  increasing  simultaneously. 
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The  vertical  velocities  are  getting  stronger  between  the  days  333  and  333.0833 
(Figure  4.18-A).  The  Langmuir  circulations  upwelling  and  downwelling  velocities  are 
very  distinguishable  in  Figure  4.18-A,  where  the  reddish  colors  represent  the  upwelling 
while  the  bluish  colors  represent  the  downwelling  regions  and  these  regions  are  detected 
in  Figure  4.18-E. 
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The  Upwelling  and  Downwelling  Points  after  Velocity  Threshold 
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Figure  4.19  Detection  of  upwelling  /  downwelling  regions  of  Langmuir  circulations 
between  the  days  333.3  and  333.3833.  The  wind  speed  is  ~1 1  ms''  and  the 
wave  height  is  ~  3  m  (the  peak  of  the  storm).  24  upwelling  regions  and  21 
downwelling  regions  (Panel  E)  were  detected. 
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The  peak  of  the  storm  event  is  shown  in  Figure  4.19  with  elean  deteeted  Langmuir 
eireulations.  The  vertieal  veloeities  in  Figure  4.19-A  are  very  strong  (about  0.06  ms"'  for 
downwelling,  0.05  ms'^  for  upwelling)  through  the  water  eolumn  at  loeations  deteeted 
Langmuir  eireulations  in  Figure  4.19-D.  The  high  wind  speed  and  energetie  waves  due  to 
this  high  wind  speed  ereates  stronger  Langmuir  eireulations  and  the  algorithm  for 
deteetion  of  these  eireulations  is  working  well. 


Figure  4.20  Summary  of  deteetion  algorithm  of  upwelling  /  downwelling  regions 
related  to  Langmuir  eireulations. 

The  diagram  in  Figure  4.20  shows  the  steps  of  the  deteetion  algorithm  of  Langmuir 
eireulations  in  a  hierarehy. 
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G.  MAXIMUM  UPWELLING  /  DOWNWELLING  VELOCITIES 

The  maximum  vertical  velocities  of  Langmuir  circulations  in  the  water  column 
were  searched  for  each  individual  detected  upwelling  and  downwelling  regions.  The 
Beam  five  vertical  velocities,  previously  low  and  high-pass  filtered  in  time,  were  4 
meters  low-pass  filtered  in  space  (vertically)  before  defining  the  maximum  velocities. 
The  depth  of  each  of  these  maximum  velocities  was  also  saved  for  the  following 
penetration  depth  analysis. 


H.  PENETRATION  DEPTHS  OF  UPWELLING  AND  DOWNWELLING 
REGIONS 

The  penetration  depths  of  each  detected  event  were  defined  as  the  maximum 
depth  where  the  coherent  Langmuir  circulations  differed  from  near-surface  values  by  a 
selected  ratio.  The  algorithm  to  define  the  penetration  depth  was  based  on  the  beam  five 
vertical  velocities  after  vertical  low-pass  filtering.  The  normalized  beam  five  velocities 
were  calculated  by  taking  the  ratio 


ratio{z,t) 


v5{z,t) 

(v5(0) 


(4.8) 


where  z  is  the  vertical  depth  and  (v5(t))  is  the  vertical  mean  of  the  upper  four  bin 

velocities  at  each  time  t.  At  each  detected  upwelling  and  downwelling  event,  this  ratio 
was  used  to  determine  the  penetration  depth  by  finding  the  depth  where  the  vertical 
velocity  had  decreased  from  the  near-surface  values  by  a  selected  ratio.  Threshold  ratios 
of  0.2,  0.3  and  0.4  were  evaluated  for  skill  in  defining  upwelling  /  downwelling  depths, 
and  three  different  bin  numbers  with  regard  to  the  three  different  thresholds  were  found 
as  representative  of  the  penetration  depth.  One  of  these  three  different  bin  numbers  was 
assigned  as  final  penetration  depth  bin  number  by  the  method  described  below.  The  bin 
number  representing  the  penetration  depth  was  converted  to  meters  since  every  bin 
number  was  0.5  m  length,  and  the  distance  between  the  reference  bin  number  and  mean 
sea  level  was  added  to  the  penetration  depth. 
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Verically  Low-pass  filtered  Beam  5  Velocities  (ms  for  Downwelling  Regions 


Vertically  Low  pas  filtered  and  Normalized  Beam  5  Velocities  for  Downwelling  regions 


Figure  4.21  Detection  of  penetration  depths  for  downwelling  regions  related  to 
Langmuir  circulations  between  the  days  352  and  352.0833  (two-hour 
period).  Panel  A  is  vertically  4  m  low-pass  filtered  beam  5  vertical 
velocities.  Panel  B  is  normalized  beam  5  vertical  velocities  by  Equation 
4.2,  the  reddish  shading  color  indicates  the  velocities  similar  to  the  near 
surface  velocities  while  the  bluish  shading  indicates  weaker  velocities  than 
the  near-surface  velocities  and  Panel  C  shows  the  penetration  depths  for 
the  downwelling  regions  after  ratio  thresholds. 


The  reddish  color  shading  in  Figure  4.2 1-B  indicates  vertical  velocities  with  values 
similar  to  near-surface  values,  while  the  bluish  color  shading  represents  the  values 
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significantly  less.  Comparisons  of  the  detected  penetration  depths  in  Figure  4.2 1-C  show 
that  the  threshold  ratio  of  0.4  eonsistently  tracks  the  depth  of  the  vertical  velocity  cells 
seen  in  Figure  4.21 -A.  0.4  threshold  ratio  was  used  to  define  the  penetration  depth  in  the 
remaining  analysis. 


Vcrically  Low-pass  filtered  Beam  S  Velocities  (m  s~ ')  for  Upwelltng  Regions 


Figure  4.22  Detection  of  penetration  depths  for  the  upwelling  regions  related  to 
Langmuir  circulations  between  the  days  352  and  352.0833  (two-hour 
period).  Panel  A  is  vertically  4  m  low-pass  filtered  beam  5  vertical 
velocities.  Panel  B  is  normalized  beam  5  vertieal  velocities  by  Equation 
4.2,  the  reddish  shading  color  indicates  the  velocities  similar  to  the  near 
surface  velocities  while  the  bluish  shading  indicates  weaker  velocities  than 
the  near-surface  velocities  and  Panel  C  shows  the  penetration  depths  for 
the  upwelling  regions  after  ratio  thresholds. 


67 


The  same  threshold  ratio  of  0.4  was  ehosen  for  finding  upwelling  cell  depth  as 
seen  in  Figure  4.22. 

I.  SPACING  BETWEEN  UPWELLING  /  DOWNWELLING  REGIONS 

Horizontal  spacing  between  successive  upwelling  regions  or  downwelling  regions 
can  be  used  to  estimate  the  horizontal  length  scale  of  the  Langmuir  circulations  described 
in  section  II. B.  For  much  of  the  time  at  the  ASIT  site,  Langmuir  circulations  were 
advected  by  strong  east- west  tidal  currents.  The  four-beam  solution  of  horizontal  velocity 
components  and  the  time  difference  between  successive  upwelling  regions  or 
downwelling  regions  were  used  to  calculate  the  spacing  during  periods  where  the 
advective  speed  exceeded  O.I  ms'\ 
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Figure  4.23  Successive  upwelling  (red  cross)  and  downwelling  (red  circles)  regions 
between  the  day  352  and  382.0833  (two  hour  period) 

For  example,  in  Figure  4.23  the  spacing  for  the  downwelling  regions  was  calculated 
between  the  consecutive  red  circles  by  multiplying  the  time  difference  between  them  with 
horizontal  advection  mean  velocity.  The  spacing  for  the  upwelling  regions  was  calculated 
between  the  consecutive  red  crosses  with  the  same  method.  This  method  requires  a  fairly 


The  Upwelling  and  Downwelling  Points  Related  to  Langmuir  curculations 
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uniform,  discreet  eell  field  to  unambiguously  define  distinet  sueeessive  upwelling  and 
downwelling  events.  This  is  generally  true  under  strong  foreing  eonditions,  but 
ehallenging  in  light  wind  /  wave  eonditions. 
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V.  ANALYSIS  AND  RESULTS 


Statistical  relationships  between  the  Langmuir  circulation  features  and  potential 
forcing  mechanisms  such  as  wind  speed  and  waves  Stokes  drift  were  explored  in  this 
long  time  series  observation.  The  tidally  forced  inner  shelf  site  provided  an  opportunity  to 
investigate  the  concurrent  effects  of  a  strong  tidal  boundary  flow  with  a  wide  range  of 
wind  forcing  conditions. 


Forcing  Mechanisms  and  Ocean  Response  with  BADCP  Data  Time  Lines 
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Figure  5.1  Wind  speed  (Panel  A),  wave  parameters  (Panel  B  and  C)  and  tidal  currents 
(Panel  D)  time  series  from  MVCO  with  the  BADCP  data  between  the  days 
228  and  360  (Panel  E)  during  CBLAST  experiment.  The  gaps  in  BADCP 
data  are  due  to  a  power  outage. 
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Winds  were  measured  4  km  north  of  ASIT  at  the  southern  shore  of  Martha’s  Vineyard, 
and  direetional  wave  speetra  1  km.  north  at  12  m  depth  to  provide  time  series  of  surface 
forcing  and  wave  condition  at  the  BADCP  site.  Southerly  winds  provide  long  fetch 
conditions,  whereas  northerly  winds  result  in  short  fetch  conditions  for  the  ASIT  location 
due  to  blocking  by  Martha’s  Vineyard  Island  seen  in  Figure  3.3.  The  tidal  currents  shown 
in  Figure  5.1-D  are  strong  for  the  experiment  region,  and  are  dominantly  in  the 
alongshore  (east-west)  direction. 

The  surface  boundary  layer  and  the  bottom  boundary  layers  are  generally 
decoupled  in  the  summer  time  due  to  strong  seasonal  stratification,  and  Langmuir 
circulations  characteristics  are  affected  by  this  stratification.  Surface  cooling  and  strong 
storm  events  diminish  the  seasonal  summertime  stratification  and  couple  the  surface  and 
bottom  boundary  layers  in  the  fall  and  wintertime.  The  available  BADCP  data  between 
days  228  and  240  are  representative  of  summertime,  and  there  are  no  significant  storm 
events  during  this  time  (Figure  5.1).  The  BADCP  data  between  the  days  280  and  360  are 
in  fall  /  wintertime  conditions  and  span  a  series  of  strong  storm  systems. 

In  this  study  locally  well-developed  seas  with  long  fetch  winds  under  the  weak 
stratification  oceanic  conditions  were  chosen  as  a  focus  for  the  analysis  of  Langmuir 
circulations.  Locally  well-developed  seas  represent  increasing  wave  conditions  with 
increasing  wind  speed  in  the  experiment  region  without  dominant  influence  of  energetic 
pre-existing  swell  from  the  Atlantic.  The  cases  meeting  these  criteria  were  selected 
between  days  280  and  360  by  plotting  the  wind  speed  time  series,  wind  /  wave  direction 
time  series,  and  time  series  of  wave  energy  density  spectrum. 
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Figure  5.2  Panel  A;  Color  plot  of  logarithmic  scale  wave  energy  density  spectrum. 

Panel  B;  Wave  and  wind  direction  time  series.  Panel  C:  Wind  speed  time 
series  between  the  year  days  284  and  304.  The  region  in  the  black  square 
is  the  selected  segment  for  this  period  of  collected  data  and  the  region  in 
the  red  circles  represent  typical  wind-induced  increasing  wave  field. 


As  the  wind  speed  in  Figure  5.2  panel  B  increases  during  year  day  288,  the  wave 
spectrum  evolves  from  short  period  local  waves  to  fully-developed  swell  and  sea 
conditions,  as  seen  in  the  red  circles  in  panel  A.  The  wind  reaches  its  maximum  speed,  18 
ms'\  on  the  day  289  with  the  most  energetic  waves  in  the  6-10  second  period  band.  After 
the  day  289,  the  winds  relax  and  the  wave  energy  decreases.  Wind  speed  is  down  to 
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2  ms''  on  day  291,  and  the  wave  energy  shows  a  trend  toward  local,  weakening, 
shortening  period  waves.  After  determining  the  locally  developed  sea  conditions,  long 
fetch  conditions  were  determined  in  the  wind  direction  time  series  (Figure  5.2-B). 
Southerly  winds  and  waves  were  classified  as  long  fetch.  The  black  square  is  the  time 
period  when  both  these  conditions  were  satisfied  and  defined  a  long-fetch,  local  wave 
field  data  segment. 
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310  315 

Time  (year  day)  -C- 


Figure  5.3  Panel  A;  Color  plot  of  logarithmic  scale  wave  energy  density  spectrum. 

Panel  B;  Wave  and  wind  direction.  Panel  C:  Wind  speed  time  series 
between  the  year  days  305  and  320.  The  Region  in  the  black  square  is  the 
selected  segment  for  this  period  of  collected  data  and  the  region  in  the  red 
circle  represents  typical  wind-induced,  increasing  wave  field  while  the 
region  in  the  purple  circle  represents  the  energetic  swell  coming  from 
Atlantic  Ocean. 
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The  red  circle  in  Figure  5.3-A  starting  on  day  317  represents  the  locally  developed  seas 
with  increasing  wind  speed.  The  region  inside  the  purple  circle  indicates  a  period  of  high- 
energy  swell  coming  from  Atlantic  Ocean.  The  region  inside  the  black  square  is  the 
selected  segment  between  the  days  314.5  and  319.  The  wind  speed  increases  up  to 
14  ms'^  during  this  time  period. 


Figure  5.4  Panel  A;  Color  plot  of  logarithmic  scale  wave  energy  density  spectrum. 

Panel  B;  Wave  and  wind  direction.  Panel  C;  Wind  speed  between  the  year 
days  324  and  342.  The  region  in  the  black  square  is  the  selected  segment 
for  this  period  of  collected  data  and  the  region  in  the  red  circle  represents 
typical  wind-induced  increasing  wave  field  while  the  region  in  the  purple 
circle  represents  the  energetic  swell  coming  from  Atlantic  Ocean. 
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The  interval  between  the  year  days  331.5  and  338.5  was  chosen  for  the  time 
period  shown  in  Figure  5.4.  Swell  from  the  Atlantic  is  very  energetic  between  the  days 
326  and  331.5  (Figure  5.4-A).  Wind  speed  starts  to  increase  to  12  ms’'  on  day  332,  and 
decreases  steadily  between  days  336.5  and  338.5. 


Figure  5.5  Panel  A:  Color  plot  of  logarithmic  scale  wave  energy  density  spectrum. 

Panel  B;  Wave  and  wind  direction.  Panel  C:  Wind  speed  between  the  year 
days  344  and  360.  The  region  in  the  black  square  is  the  selected  segment 
for  this  period  of  collected  data  and  the  region  in  the  red  circle  represents  a 
typical  wind-induced  increasing  wave  field.  The  purple  square  region  is  a 
one-day  period  between  the  days  351  and  352  for  stratification  and  shear 
profile  analysis. 
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There  are  four  clear  long  fetch  storm  cases  in  the  day  344  to  360  interval  seen  in  Figure 
5.5,  so  the  whole  time  period  shown  in  Figure  5.5  was  chosen  for  the  long-fetch,  low- 
stratification  analysis.  The  purple  region  represents  a  one-day  period  when  the  wind 
speed  increased  from  1  ms’'  to  12  ms'^  between  the  days  351  and  352,  and  the  current 
profile  was  examined  as  follows  for  this  period. 
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Figure  5.6  Panel  A;  Vertically  averaged  barotropic  alongshore  (u)  and  cross-shore 
velocities  (v).  Panel  B;  Demeaned  alongshore  velocity  profile  through  the 
water  column.  Panel  C;  Demeaned  cross-shore  velocity  profile  through  the 
water  column  between  the  year  days  351  and  352.  Vertical  blue  line 
represents  the  condition  of  low  wind  and  high  tidal  current  during  351.1 
and  351.14,  vertical  red  line  represents  the  condition  of  low  wind  and  low 
tidal  current  during  351.32  and  351.36,  vertical  blue  line  represents  the 
condition  of  high  wind  and  tidal  current  during  351.96  and  352. 


Vertical  mean  horizontal  velocity  time  series  and  demeaned  current  profile  time 
series  are  shown  in  Figure  5.6  with  along-shore  (east-west)  tidal  currents  greater  than 
cross-shore  tidal  currents  at  this  site  for  a  24  hour  period.  During  this  time  wind  speed 
before  the  351.6  is  between  0  and  2  ms'*,  then  increases  to  12  ms'*  on  day  352.  Velocity 
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differences  between  the  upper  and  bottom  ocean  are  clearly  seen  in  the  demeaned 
velocity  profdes  in  Figure  5.6-B  and  Figure  5.6-C,  under  low  wind  conditions  (before  day 
316.5).  Once  the  wind  speed  starts  to  increase  after  day  351.5,  the  pre-existing  velocity 
shear  rapidly  weakens  through  the  water  column.  This  suggests  that  the  wind 
homogenizes  existing  shear  under  these  conditions.  During  cold  weather,  strong  wind 
forcing  conditions  in  the  fall  at  this  site,  the  water  column  is  largely  unstratified.  During 
these  conditions,  the  strong  tidal  currents  are  likely  to  form  a  Taw  of  the  wall’ 
logarithmic  velocity  profile  with  constant  stress  across  the  water  column.  This  hypothesis 
was  tested  by  plotting  a  logarithmic  depth  scale  versus  horizontal  velocities  in  the  water 
column  for  a  combination  of  high  and  low  winds  and  tidal  currents. 
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Figure  5.7  Logarithmic  Depth  scale  and  horizontal  velocity  plots  to  obtain  the 
stratification  for  three  different  conditions  described  in  Figure  5.6.  The 
colors  are  the  same  with  the  colors  used  in  Figure  5.6  to  represent  the 
three  different  conditions. 

The  linear  decrease  of  horizontal  velocities  with  respect  to  logarithmic  scale  indicates  a 
well-developed  constant  stress  flow  at  low  wind  speeds.  The  dashed  red  and  blue  lines  in 
Figure  5.7  represent  the  low  wind  speed  conditions,  while  the  black  dotted  line  shows 
high  wind  condition,  using  the  same  colors  in  Figure  5.6.  The  low  wind  and  low  current 


Logarithmic  Depth  vs  Honzontal  Veloctty 
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condition  profile  in  Figure  5.7  decreases  almost  linearly  through  the  water  eolumn, 
indieating  an  unstratified  weaker  tidal  boundary  layer.  After  wind  speed  inereases  to  12 
ms'*  the  dotted  black  line  in  Figure  5.7  shows  this  log-layer  to  have  been  replaeed  with  a 
low  shear,  homogenized  water  eolumn. 
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Figure  5.8  Demeaned  Horizontal  velocities  vs.  Depth.  Panel  A;  The  demeaned 
horizontal  veloeity  profile  through  the  water  column  for  high  tidal  currents 
(greater  than  0.1  ms'*).  Panel  B;  The  demeaned  horizontal  velocity  profile 
for  low  tidal  currents  (less  than  0.1  ms'*). 

These  profiles  are  shown  as  demeaned  linear  z  profiles  in  Figure  5.8  to  emphasize  the 
differences  in  shear  structure  as  wind  increases.  In  Figure  5.8-A  with  high  tidal  eurrent, 
the  wind  speed  for  the  blue  dashed  line  shear  profile  is  about  2  ms'*  between  the  days 
351.1  and  351.14  while  the  wind  speed  for  red  erossed  line  shear  profde  is  about  8  ms'* 
between  the  days  351.55  and  351.59.  Demeaned  horizontal  velocity  shear  with  2  ms'* 
wind  speed  is  greater  than  demeaned  horizontal  velocity  shear  with  8  ms'*  wind  speed, 
while  the  dotted  blaek  line  shear  profile  with  12  ms'*  wind  shows  very  low  velocity  shear 
for  the  time  between  the  days  351.96  and  352.  In  Figure  5.8-B  with  low  tidal  currents,  the 
blue  dashed  line  demeaned  horizontal  velocity  shear  profile  between  the  days  351.32  and 
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351.36  with  2  ms"'  wind  speed  is  much  greater  than  the  dotted  black  line  demeaned 
horizontal  velocity  shear  profile  between  the  days  351.82  and  351.86  with  10  ms"'  wind 
speed.  For  both  high  tidal  and  low  tidal  current  cases,  the  bottom  boundary  layer  shear 
profiles  in  the  water  column  rapidly  decrease  with  increasing  wind  speed  while  surface 
forced  mixing  processes  increase. 


CBLA3T  Wind  Speed,  SWH,  Up/Dov/nwelling  Velocities 


Figure  5.9  Panel  A;  Time  series  of  wind  speed.  Panel  B;  Significant  wave  height 
(SWH)  time  series.  Panel  C;  Detected  downwelling  /  upwelling  regions 
related  to  Langmuir  circulations  responsible  for  the  mixing  process 
between  the  days  351  and  352 
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A  one-day  time  series  of  wind  speed  (Figure  5.9-A)  and  signifieant  wave  heights  (Figure 
5.9-B)  show  inereasing  wave  heights  with  the  inereasing  wind  speed.  The  low  wind 
speed  and  low  signifieant  wave  heights  before  the  day  351.5,  result  in  low  levels  of  wind 
and  wave  interaetion,  henee  the  amount  of  the  deteeted  downwelling  /  upwelling 
veloeities  related  to  the  Langmuir  eireulations  is  small.  After  day  351.5,  the  wind  and 
wave  interaetion  effeet  is  inereasing,  and  ereating  more  Langmuir  eireulations  as  shown 
in  Figure  5.9-C  above  our  deteetion  thresholds.  We  hypothesize  that  the  homogenizing 
proeess  aeross  the  water  eolumn  after  351.5  is  due  to  the  inereasing  strength  of  the 
Langmuir  eireulations  aeting  in  eoneert  with  surfaee  shear-driven  turbulenee. 

The  days  between  288  and  295.5;  314.5  and  319;  331.5  and  338.5;  344  and  360 
were  seleeted  to  represent  unstratified,  loeally  well-developed  seas  with  southerly,  long 
feteh  wave  /  wind  eonditions.  It  is  expeeted  that  strong  Langmuir  turbulenee  oeeurs  for 
these  time  periods.  Langmuir  number  is  a  parameter  to  prediet  if  the  Langmuir  turbulenee 
is  strong  or  not.  The  Langmuir  number  was  first  deseribed  by  Leibovieh  (1977a)  as 
shown  in  Equation  2.38.  MeWilliams  et  al  (1997)  parameterized  the  turbulent  Langmuir 
number  as, 

La,={uJUj'^  (5.1) 

where  u,  is  the  water  frietion  veloeity  and  the  is  the  surfaee  Stokes  drift  veloeity. 

Water  frietion  veloeity  is  parameterized  through  a  widely  used  empirieal  drag  eoeffieient 
by  Large  and  Pond  (1981)  as 

M,  =  (r/  yC>  )*^^ 

^  =  PairCDUm  (5-2) 

=(4.9x10^ +6.5x10^' xMio) 

where  r  is  the  wind  stress  at  the  sea  surfaee,  is  the  drag  eoeffieient,  Mjq  is  the  wind 
speed  measured  at  10  m  above  the  sea  surfaee  and  p  is  the  density.  The  surfaee  Stokes 
drift  is  defined  by  Philips  (1977)  in  Equation  5.3. 
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(5.3) 


=  (jka^ 

CT  =  [gkf^ 

where  k  is  the  dominant  swell  wave  number,  a  is  the  wave  amplitude,  g  is  the 
gravitational  acceleration  and  cris  the  particular  wave  frequency  associated  with 
dominant  swell  wave  number. 

If  the  surface  Stokes  drift  velocity  is  low  for  a  given  wind  forcing,  the  turbulent 
Langmuir  number  will  be  large  and  shear  turbulence  will  dominate.  The  smaller  the 
turbulent  Langmuir  number  (stronger  Ug),  the  stronger  the  Langmuir  circulation.  Li  et  al 
(2004)  defined  a  turbulent  Langmuir  number  of  approximately  0.7  for  the  transition  from 
shear  turbulence  to  Langmuir  turbulence. 

Downwelling  Velocities  vs.  Langmuir  Number 
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Figure  5.10  Downwelling  velocities  and  Langmuir  number  relationship  between  the 
days  288  and  295.5;  314.5  and  319;  331.5  and  338.5;  344  and  360. 
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Upwclling  Vclocitic$  vs.  Langmuir  Number 
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Figure  5.11  Upwelling  velocities  and  Langmuir  number  relationship  between  the  days 
288  and  295.5;  314.5  and  319;  331.5  and  338.5;  344  and  360. 

Both  Figure  5.10  and  5.11  show  a  strong  increase  of  detected  downwelling  /  upwelling 
velocities  largely  occurring  with  a  range  of  Langmuir  numbers  from  0.1  to  0.4,  which 
confirms  that  the  Langmuir  turbulence  is  strongest  for  low  Langmuir  numbers  where  the 
wave  Stokes  drift  is  high.  The  velocities  shown  with  the  Langmuir  number  greater  than 
0.5,  largely  represent  short  fetch  within  the  “long  fetch”  time  series.  In  these  cases, 
surface  shear  stress  dominates  over  Langmuir  turbulence.  The  period  between  days  291 
and  293  is  very  good  example  for  these  cases  (Figure  5.12  and  Figure  5.13).  In  Figure 
5.2-C,  winds  are  northerly  and  wind  speed  is  ~2  ms'^  on  day  291  and  increases  up  to  7 
ms'*,  then  decreases  to  1  ms'*  between  the  days  291  and  292,  then  starts  to  increase  and 
reaches  6  ms'*  on  day  293.  Under  short  fetch  conditions  the  Langmuir  number  is  greater, 
and  the  vertical  velocities  are  lower  for  this  time  period.  Li  et  al  (2004)  suggested  a 
Langmuir  number  of  10,  which  results  in  a  strong  shear  turbulence  and  very  small 
Langmuir  turbulence,  but  this  is  not  seen  for  the  forcing  conditions  seen  here. 
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Figure  5.12 


Figure  5.13 


Downv/elling  Velocited  vs.  Langmuir  Number  between  days  291  and  293 


Downwelling  veloeities  with  Langmuir  number  between  the  days  291  and 
293  with  northerly  winds  (short  fetch) 
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Downwelling  velocities  with  Langmuir  number  between  the  days  291  and 
293  with  northerly  winds  (short  fetch) 
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Figure  5.14 


Figure  5.15 
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Downwelling  velocities  and  friction  velocity  regression  plot  for  tidal 
currents  greater  than  0.1  ms'^  for  the  selected  days. 


Downwelling  Velocities  vs.  u.  for  Low  Tidal  Current 
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Downwelling  velocities  and  friction  velocity  regression  plot  for  tidal 
currents  less  than  0.1  ms"'  for  the  selected  days. 
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Figure  5.16  Upwelling  velocities  and  friction  velocity  regression  plot  for  tidal  currents 
greater  than  0.1  ms''  for  the  selected  days. 
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Upwelling  Velocities  vs  u.  for  Low  Tidal  Current 
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Figure  5.17  Upwelling  velocities  and  friction  velocity  regression  plot  for  tidal  currents 
less  than  0.1  ms"'  for  the  selected  days. 

Statistical  relationships  between  the  water  friction  velocity  and  downwelling  / 

upwelling  velocities  for  high  tidal  currents  and  low  tidal  current  were  examined  for  the 
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selected  days  to  see  whether  the  presence  of  a  strong  tidal  boundary  layer  through  the 
water  column  affected  generation  of  Langmuir  cells.  Correlations  between  downwelling 
velocities  and  surface  friction  velocity  are  statistically  significant  for  both  high  tidal 
currents  and  low  tidal  currents.  In  Figure  5.14  for  high  tidal  current,  the  downwelling 
velocities  increase  with  increasing  friction  velocity  and  the  regression  coefficient  is  0.44 
with  a  slope  of  1.05.  In  Figure  5.15  for  low  tidal  current,  the  regression  coefficient  is  0.41 
with  a  slope  of  1 .09  between  the  downwelling  velocities  and  friction  velocity,  suggesting 
no  differences  due  to  the  presence  of  a  strong  tidal  bottom  boundary  layer. 

Similarly,  upwelling  velocities  and  friction  velocities  are  statistically  significant 
for  both  high  tidal  currents  and  low  tidal  currents.  In  Figure  5.16  for  high  tidal  current, 
the  upwelling  velocities  are  increasing  with  increasing  friction  velocity  and  the  regression 
coefficient  is  0.50  with  a  slope  of  1.1.  In  Figure  5.17  for  low  tidal  current,  the  regression 
coefficient  is  0.45  with  a  slope  of  1.1  between  the  upwelling  velocities  and  friction 
velocity. 

The  linear  relationship  found  between  the  friction  velocity  and  downwelling 
velocities  is  consistent  with  results  of  Leibovich  (1983)  shown  in  Figure  2.10.  We  also 
found  a  similar  linear  relationship  for  upwelling  velocities  in  addition  to  his  result.  The 
slope  of  the  regression  line  is  -1.07  for  downwelling  velocities  while  it  is  -1.1  for 
upwelling  velocities,  which  indicates  that  there  is  a  slight  asymmetry  between  them. 
Since  wind  and  wind  generated  waves  are  responsible  for  the  formation  of  these 
circulations,  the  strength  of  the  circulations  should  increase  with  increasing  forcing 
strength.  The  regression  coefficients  are  not  high,  but  we  hypothesize  that  this  is  because 
of  the  wide  range  of  scales  of  the  Langmuir  circulations  within  the  water  column  for 
moderate  and  strong  forcing.  Langmuir  (1938)  and  many  others  also  observed  that  small- 
scale  circulations  exist  within  large-scale  circulations. 
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Downwelling  Penetration  Depth  vs.  u.  for  High  Tidal  Current 


Figure  5.18  Downwelling  Penetration  Depths  and  friction  velocity  regression  plot  for 
tidal  currents  greater  than  0.1  ms"'  for  the  selected  days. 

The  relationship  between  downwelling  penetration  depths  and  the  friction 
velocity  is  shown  in  Figure  5.16  with  a  regression  coefficient  of  0.19  and  they  are  not 
statistically  correlated.  This  confirms  that  penetration  depth  of  the  Langmuir  circulations 
does  not  depend  on  the  forcing  mechanism  but  rather  is  related  to  water  depth  and 
stratification  profiles  at  the  site. 
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Figure  5.19  Downwelling  penetration  depth  and  spacing  (between  the  downwelling 
regions)  relationship  for  days  348  and  349;  351  and  352;  355  and  356, 
which  represents  increasing  wind  and  wave  forcing. 
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The  correlation  between  cell  spacing  and  penetration  depth  was  found  for  both 
increasing  wind  speed  and  wave  heights,  and  constant  high  winds  and  wave  heights.  The 
regression  coefficient  of  0.44  with  a  slope  of  8.7  was  found  for  increasing  wave  height 
and  wind,  which  indicates  a  statistical  correlation  (Figure  5.19).  Nevertheless,  the  same 
correlation  was  not  found  for  the  high  constant  wind  speed  and  wave  conditions  which 
are  representative  of  the  older  waves  case  with  a  small  Langmuir  number  of  ~  0.3. 
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Figure  5.20  Downwelling  penetration  depth  and  spacing  (between  the  downwelling 
region)  relationship  for  the  days  selected  for  the  days  349  and  350;  352 
and  353;  356  and  357,  which  represents  high  wind  speed  and  waves. 

The  regression  coefficient  of  0.15  in  Figure  5.20  with  slope  of  4  indicates  that 
there  is  not  a  significant  correlation  between  spacing  and  downwelling  penetration  depths 
for  high  wind  and  wave  conditions,  again  suggesting  a  wide  range  of  cell  scales  for  these 
conditions. 

The  detection  of  Langmuir  circulations  under  stratified  summer  ocean  conditions 
is  complicated  by  the  presence  of  strong  internal  waves.  Fligh  vertical  velocity  internal 
solitons  might  have  coherent  vertical  velocities  like  Langmuir  circulations  in  the  upper 
ocean.  Evidence  of  this  can  be  seen  in  the  summer  time  velocity  profiles  between  days 
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228  and  242.  Strong  velocity  shear  can  be  seen  in  a  24-hour  demeaned  current  profile 
time  series,  shown  in  Figure  5.21  panel  B  and  panel  C. 


Barotropic  Horizontal  Velocities 

0.8  r  I  I  I - 1  1 - 1  I - r 


Demeaned  Along-shore  (U)  velocities  (m  s~ ') 


Demeaned  Cross-shore  (V)  velocities  (m  s~^) 


231.5  231.6  231.7  231.8  231.3  232  232.1  232.2  232.3 

Time  (Year  Day) 

-C- 


Figure  5.21  Panel  A;  Vertically  averaged  barotropic  along-shore  (u)  and  cross-shore 
velocities  (v).  Panel  B;  Demeaned  along  shore  velocity  profile  through  the 
water  column.  Panel  C:  Demeaned  cross-shore  velocity  profile  through  the 
water  column  between  the  year  days  231.4  and  232.4  during  the  summer 
time  under  strong  stratification  condition. 

The  wind  speed  is  up  to  8  ms'\  but  mixing  is  constrained  to  the  upper  few  meters  by 
strong  stratification  during  this  time  period.  The  internal  solitons  having  strong  and 
coherent  velocities  are  shown  in  Figure  5.22  by  black  circles.  These  were  detected  as 
Langmuir  circulations  (Figure  5.22-E).  The  analysis  of  Langmuir  circulations  detected 
under  these  conditions  may  give  inaccurate  results.  For  example,  we  found  penetration 
depths  down  to  12  meters  for  detected  upwelling  /  downwelling  events  in  Figure  5.22-E, 
which  are  unlikely  for  16.5  m  depth  of  water  with  strong  summer  stratification.  Erom 
Eigure  5.21-B,  the  penetration  depths  of  the  Langmuir  circulations  should  be  a  maximum 
of  3  to  5  m.  Shallow  Langmuir  circulations  constrained  by  near-surface  stratification  are 
not  well  resolved  by  the  five  beam  BADCP  and  analysis  techniques  used  in  this  study 
due  to  the  co-existing  strong  internal  solitons  during  the  summer  time. 
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Figure  5.22  Panel  A:  Vertical  velocities  trough  the  water  column.  Panel  B;  Upper  four 
bin  vertical  velocities.  Panel  C;  Vertical  mean  of  the  velocities  in  panel  B 
and  all  the  upwelling  (+)  /  downwelling  (-)  events.  Panel  D;  Remaining 
upwelling  /  downwelling  events  after  velocity  thresholds.  Panel  D; 
Detected  upwelling  /  downwelling  events  as  Langmuir  circulations 
between  the  days  232  and  232.0833.  The  black  circles  in  Panel  A  are 
internal  solitons,  which  are  detected  as  Langmuir  circulations. 
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VI.  CONCLUSIONS 


Analyses  of  detailed  water  eolumn  veloeity  profiles  for  Langmuir  eireulations  at 
an  inner  shelf  site  foeused  on  long-feteh,  unstratified  water  eolumn  with  loeally 
developed  sea  eonditions.  A  deteetion  algorithm  for  upwelling  /  downwelling  veloeities 
due  to  Langmuir  eireulations  was  ereated  using  high-resolution  vertieal  veloeity  profdes 
from  a  five  beam  BADCP  deployed  during  CBLAST-LOW  at  the  ASIT  researeh  tower 
south  of  Martha’s  Vineyard. 

The  observations  show  that  a  wind-foreed  surfaee  layer  that  generates  strong 
Langmuir  eireulations  greatly  reduees  water  eolumn  shear  arising  from  a  tidally  foreed 
bottom  boundary  layer.  Langmuir  numbers  were  found  between  0.1-1. 5  for  the  time 
period  of  deteeted  upwelling  /  downwelling  veloeities  of  Langmuir  eireulations,  but  the 
strongest  vertieal  veloeities  were  found  with  the  Langmuir  number  of  0.1  to  0.4.  This  is 
eonsistent  with  the  modeling  results  of  Me  Williams  et  al  (1997),  who  found  Langmuir 
number  0.3  most  favorable  for  Langmuir  eireulations  to  develop.  The  short  feteh 
eonditions  with  Langmuir  number  greater  than  0.5  have  weaker  vertieal  veloeities  than 
the  long  feteh  eonditions  with  Langmuir  number  less  than  0.5.  The  Langmuir  eireulation 
effeets  on  the  mixing  proeesses  inerease  with  inereasing  wind  speed  and  wave  heights. 

The  number  and  strength  of  Langmuir  eireulations  inerease  with  inereasing  wind 
speed  and  wave  heights.  Downwelling  veloeities  were  found  from  0.01  to  0.08  ms'^ 
while  upwelling  veloeities  were  found  from  0.05  to  0.07  ms'^  with  a  small  differenee, 
whieh  indieates  a  small  asymmetry  between  them. 

The  regression  eoeffioients  found  between  water  frietion  veloeity  and  both 
downwelling  /  upwelling  veloeities  show  that  Langmuir  eireulation  strength  is 
statistieally  eorrelated  with  the  joint  surfaee  stress  and  Stokes  drift  foreing  meehanisms. 
The  same  eorrelations  were  found  with  high  tidal  and  low  tidal  eurrent  eonditions, 
suggesting  that  there  is  no  signifioant  effeet  from  the  strong  tidal  boundary  layer  on  the 
formation  of  Langmuir  eells.  Co-existenee  of  multiple  seale  Langmuir  eireulations  results 
in  low  regression  eoeffioients  (i.e.  from  0.4-0. 5),  but  the  regressions  still  represent  a 
statistieally  signifioant  oorrelation. 
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The  penetration  depths  of  deteeted  eireulations  were  found  to  be  from  6  to  14 
meters  and  independent  of  the  foreing  strength.  Langmuir  spaeing  between  the  eells  was 
found  to  be  from  10  to  160  m,  and  a  statistieally  signifieant  regression  eoefficient  of  0.44 
was  found  between  the  spaeing  and  penetration  depths  for  inereasing  wind  speed  and 
wave  heights.  In  eontrast,  the  statistieal  eorrelation  deereased  to  0.15  for  eonstant  high 
wind  and  wave  eonditions  with  small  Langmuir  number. 

Langmuir  eells  are  energetie,  vertieally  eoherent,  large  horizontal  seale 
eireulations  within  the  oeean  surfaee  layer  that  aet  to  eliminate  vertieal  gradients  of 
eurrent  and  other  properties.  Improved  parameterizations  of  Langmuir  eireulation  are 
being  investigated  in  eollaboration  with  ongoing  LES  modeling  by  Jim  MeWilliams  and 
Ming  Su. 
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